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ABSTRACT 

A computer program that can determine the collapse load 
of a transversely loaded grillage that allows for various 
types of loading and failure modes is presented- The 
program allows the user to constrain particular grillage 
members to have fully developed plastic hinges at the 
location of maximum moment which gives the program the 
capability to study redundancy of transversely loaded 
grillages- The lower bound theorem of limit analysis and 
linear programming techniques are used to determine a grid 
load factor and provide a failure mode and associated 
bending moments- A user's guide and program explanation are 
provided- An example grillage is studied using this program 
and the load factor and redundancy is discussed- 
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INTRQDUCTIQN 

The capabilility to easily and effectively study a 
structure from the stand point of redundancy would 
significantly contribute to the tools available to engineers 
for calculating safety factors and provide regulatory bodies 
an inexpensive way to verify or develop safety factors with 
analytical basis. 

Strengthening of ship structures has been an area of a 
significant amount ot work- Ice loading is particularly 
threatening to ship safety and also has the character i sti c 
of very large localized or point loads- The program 
presented here is ideal for the study of local loading or 
fail ure- 

Previous work done on grillages that used similar 
techniques to determine the load carrying capability and 
possible failure modes was done by Tait and Hodge (1) in 
their development of a computer program called BEAMPLT that 
utilizes linear programing techniques to find a solution to 
transversely loaded grillages. This work was later expanded 
upon to provide a simple method of computing failure modes 
and load factors under various patch loads by Abbott (2) - 

Hodge (1), in his paper on transversely loaded 
grillages, demonstrated the feasibility of using linear 
programing to calculate a failure mode and load carrying 



“ 8 - 



capacity (load factor) for transversely loaded grillages- 
His work allowed for several types of boundary conditions. 
These conditions included simply supported, clamped, and 
cantilevered. The loading scheme was limited to point loads 
on each node. 

Abbott (2) used this approach to solve for various 
types of patch loads on a grillage. He wrote a computer 
program that employed the same techniques as Hodge used. 

He employed linear programing techniques for calculating a 
failure mode and load factor for a given grillage. The 
patch loads that he was interested in were those that could 
be associated with relatively small area ice loading on the 
side of a ship. These patch loads were transferred to the 
nodes of the grillage using a simple lever principle. 

Both Hodge (1) and Abbott <2) allowed failures to occur 
only at nodes in the grillage. Hodge’s scheme allowed 
loading only at the nodes. Abbott's method was adequate to 
estimate failure with reasonable accuracy, but his 
assumptions in transferring patch loads into nodal loads was 
simplistic and did not allow for point loads off of a node. 
Abbott's approach followed Hodges very closely and much 
programing effort was made to allow for all the boundary 
conditions that Hodge had programmed for- 

The computer program presented here is simplified for 
implementation on a VAX computer with one of the standard 
Internat 1 onal Mathematical and Statistical Library's linear 
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programing routines- Linear programming is used because it 
is fast, inexpensive and easily adaptable to many 
geometries. I have allow for only those boundary conditions 
that are frequently encountered in ship structures, that is 
clamped or simple supported- The capability to handle point 
or linear line loads on any grid member has been introduced. 
The ability to handle a uniform load over the entire 
grillage is also important and should be developed- The 
most important aspect , however, is the capability of adding 
additional constraints to the linear programing routines 
that will allow for failures at other than a node. An 
iterative approach to obtain a final solution to this 
general loading scheme is employed. One of the many types 
of problems that this approach is well suited for is 
redundancy study- In developing the capability to introduce 
local failures and the ability to handle large local loads 
the program finds particular suitable use in analyzing ice 
loaded ships structures or slow collisions- I will 
demonstrate its applicability for this by discussing the 
redundancy of a grillage of a ship's side that is loaded as 
one might expect to encounter in ice loaded ship structures. 
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CHAPTER 1 

MODELING CONSIDERATIONS 

The problem of transversely loaded grids covered with 
plating is a complex one that, except for very simple cases 
does not have a closed solution. There are several aspects 
of this problem that can be easily simplified without 
significantly effecting the end results. In ship structures 
it is the goal to obtain a reasonable bound for the ultimate 
strength of a given structure in order to minimize its 
weight and in doing so maximize its efficiency. It is well 
known that the thin skin plating adds little to the 
transverse strength of a grillage and can easily be included 
as an effective width calculation in determining the fully 
plastic bending moment of the individual members of the 
grillage. As such, this paper will deal entirely with the 
strength of the grid. When a grid is deflected, twisting 
moments are introduced in grid members. For small 
def or mati ons, these twisting deformations are relatively 
small and will be ignored. It is also assumed that shear 
force does not contribute to the breakdown at any plastic 
hinge. For analysis, all grid members are assumed to be 
perfectly elastic and undergoing only small deformations. 

The program will identify a lower bound for ultimate 
collapse and give a possible failure mode. 

The problem reduces to a limit analysis of grillages 
loaded with a variety of transverse loads- In order to 
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apply linear programming techniques each load must be 
reduced to correspond! ng point loads at nodes. This 
reduction does not adversely affect the outcome since it is 
at these intersections that the interactions occur. The 
behavior of the member itself is determined later using the 
moments calculated at each node as the end moments for the 
given member. Once the load reduction is complete the lower 
bound plastic load limit is found using a standard linear 
programmi ng rout i ne. 

As a simple example consider a simple supported nine 
node grid consisting of six members arranged sysmmetr i cal 1 y 
as shown in figure (1). Applying a concentrated load F at 
node number 5, the center node, will produce a failure mode 
as shown. The circles indicate fully developed plastic 
hinges. Twisting does not occur in this example because of 
symmetry. This is not the general case. The effect of 
torque in normally encountered structures is relatively 
small when compared to bending, and in ignoring twisting 
moments results will be conservative. Tait <1) has shown 
that for typical I sections the error is less than 0.1 
percent for the types of grids generally considered- For 
long narrow less symmetrical grids the error may be as large 
as S percent - 

Determining the work done in creating the deformation 
pattern shown in figure 2., which is the only symmetrical 



failure mode possible for this geometry, the lower plastic 



N 



- 12 - 




Figure 1- Typical grillage with general loading 




1 




Figure 2. A nine node grid collapse 
and deformation profiles 



mode 
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1 i fni t for collapse can be determined- If the angle of 
rotation is G then the transverse load is displaced by 0L 
where L is the length of each member. And, the external 
work is F0L. Each hinge shown undergoes a rotation of 9 so 
the total work performed is 4Mo9. By setting the two values 
equal to each other the total collapse load can be 
cal cul ated . 



FL9 = 4Mo9 
F = 4Mo/L 

The program GRIDS presented here produces identical 
solutions- The results of other more complex examples 
completed by Hodge (3) and Tait (1) have also been 
duplicated and the results have been identical- GRIDS is a 
FORTRAN computer program that uses plastic methods of 
analysis to determine the response of a grillage subjected 
to a general scheme of transverse loading- The program can 
handle linear line loads on a member, point loads at any 
location on the grillage or a distributed load over the 
entire grillage- This program can be used to find the 
plastic load limit of any grillage- 

GRIDS is a relatively simple program that employs a 
linear programing technique to solve a system of equations 
that are constrained by the limit of a fully developed 
plastic moment in a member- A program developed by Tait and 



- 14 - 



Hodg 0 (1) producBd a similar solution but was limited in the 
type of loading that it was able to handle. It allowed only 
point loads at beam intersections- Abbott (2) modified 
Tait's and Hodge's program and introduced a simple geometric 
method of reducing patch loads to nodal loads (loads at 
intersections) and then applied Tait's and Hodge's method 
directly to arrive at a solution. GRIDS uses the same 
technique combined with the specific loading scheme to check 
individual grillage members for potential local failures and 
then employ an iterative process in narrowing the bound and 
failure mode based on the points in the grillage not on the 
nodes that are most likely to fail. 

Program Input and Problem Approach 

The input data used by GRIDS consists of the specific 
geometry, number of beam intersections, beam yield moments, 
lengths of the beams, boundary conditions, and specific data 
for all the loads. The program is limited to one point load 
per node, one additional point load per member and one 
linear line load per member. 

Once all the loading data is entered, equations of 
equilibrium are developed for each member. In reducing 
these equilibrium equations for each member the assumption 
is made that no local failure of the member has occurs- 
This is verified with the output from the first iteration of 
GRIDS. If a particular member is found that has exceeded 
its yield moment then the moment at the location of this 
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Figure 



Free Body Diagram for a Node and Surrounding 
members (general loading condition) 
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occurrence IS constrained. An additional equation is added 
to the system for each such member. A new solution is 
determined which includes these new constraints. Figure (3) 
shows an equilibrium diagram for a node in a general loading 
condi t i on . 

In order to calculate a solution the following approach 
is followed- The shear forces at the ends of each member is 
solved for and a moment balance performed. For example, 
taking member (i) we have: 
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The total shear at a node will be the sum of the shears of 
the members that meet at that node- The only unknowns in 
these equations are the end moments. Combining the load 
terms and end shears into a single load for each node f we 
solve for the vertical and horizontal moments at each node- 
Ref erring to figure (3) we have from equilibrium of 
force at node ij ; 



f . = V2 + V3 + V6 + V7 

IJ 



(3) 



and from equilibrium of moments across each member we have 



V3 = V4 = i>r (4) 

L 

Similar expressions are obtained for all other members. 
Substituting equation <4) into equation (3) we have 



/*i j-*^i jH-1 H 

S-1 



li. . M. 

L 

J 

PP. f . . 
1 J 1 J 



M. M , . 
11 - 1 -1 1 

4-1 



(5) 



which satisfies equilibrium of node ij- Similar equations 
are developed for each node. 

The program GRIDS constrains each moment to satisfy 



-18- 



-Mo < M. . < Mo (6) 

1 J 

-Mo' < M' . . < Mo' <7) 

1 J 

where Mo is the yield moment of a particular member- The 

yield moments of symmetrical beams can easily be calculated 

from the geometry and material properties of the beam- 

Now GRIDS applies a linear programming routine to the 

system to find a solution. To do this the program will 

maximize the load factor PP, where PP is the multiplier of 

each of the applied loads required to induce collapse. In 

maximizing PP a maximum number of moments M. and M' will 

1 J 1 J 

reach the imposed constraint j: Mo- 

The linear programming problem seeks to minimize an 

✓ 

objective function that is a linear function of unknowns, 
subject to constraints. These constraints consist of linear 



equal i t i es 


and inequalities. In 


the standard form the 


inequalities can be 


expressed as 


equal i ties 


and 


the problem 


becomes: 
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The coefficients of the constraint equations and of the 
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objective function are known values and the X.'s are 

1 



unknowns. 

Since our problem is to maximize PP it should be noted 
that minimizing f (XI , X2, X2, . . . . , Xn) is equivalent to 
maximizing -f ( X 1 , X2, X3, . - - - ? Xn ) - A detailed description of 
linear programing can be found in reference (4). 

The objective function of our problem becomes: 




+ 



C X = PP 



( 10 ) 



n n 



where the X's are given by 



X 



M. + M 
1 j 



( 11 ) 



n 



o 



so that the inequalities (6) and (7) become 



0 < X < 2M 



( 12 ) 



n o 



o 



O < < 2M^a 



(13) 
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CHAPTER 2 
PROGRAM OPERATION 

GRIDS is a Fortran program written in a modular fashion 
that makes maximum use of separate subroutines and common 
memory. This approach makes the program much easier to 
understand and gives it the flexibility to be easily 
changed . 

The program has some specific limitations- All 
horizontal members and their maximum yield moments must be 
the same. All vertical members and their maximum yield 
moments must be the same. However, vertical and horizontal 
members need not be the same- Boundary conditions may only 
be clamped or simple supported- But, the boundary condition 
of each side is independent of the others- Point loads may 
be applied at any point of on the grillage, but there must 
be only one point load per member with one additional point 
load allowed at each node. One linear line load per member 
is allowed. All loads must be applied perpendi cul ar and 
transverse to the grillage. 

The main program GRIDS assimilates major variables and 
input and sequences the subroutine operation. There is only 
one set of significant calculations performed in the main 
program. The calculation consists of determining the 
dimensions of the arrays that the IMSL linear programming 
routine ZX3LP uses and could not be done in the subroutine 



that actually calls 2X3LP. 
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Each of the program's functions is controlled by an 
independent subroutine. The controlling program GRIDS takes 
the user input for the grillage to be studied and assigns 
variables for the boundary conditions and geometry of the 
structure. Subroutine INTPL initializes the load array. 

MINT uses the specific boundary conditions and initializes 
the output load arrays HMNT ( i ) and VMNT(i). The index i is 
for the i th node. These arrays are filled with ones and 
zeros. All nodes that are known to have a zero moment 
because of the boundary conditions are assigned a zero, all 
other unknown moments are assigned a one- These arrays are 
used to determine the total number of unknowns and are used 
as multipliers in Subroutine CONST to ensure that the 
boundary conditions are handled properly (ie- zero moments 
at simple supports). Later they are used to store the nodal 
moments calculated in the linear programming routine- 

Subroutine ENTRL handles assigning all loads to the 
proper arrays and decomposes them into nodal loads for use. 
Each load must be stored so that individual beam moment 
distributions can later be calculated by Subroutine BEAMCK. 

Once the loads are entered the user decides if yield 
moments are allowed to form in the members away from the 
nodes- If he chooses to allow additional plastic hinges 
then the subroutine CONST is called and the additional 
constraint equations are written. These equations may be 
generated automati cal 1 y in which case they are based on the 
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grillage scheme and an initial check of each member. They 
can also be generated one at a time manually. In this case 
the user inputs which member is to have the additional 
constraint and the location on the member that is expected 
to develop the plastic hinge. In addition, if the manual 
mode IS chosen, the constraint equation can be written 
normally as a likely failure point or it can be written as 
if the point has already failed due to some very large load 
(in this case it becomes an equality constraint). These 
possibilities are the heart of the program and are discussed 
in detail in the body of this paper. 

Subroutine EDNS is used to calculate the nodal 
equilibrium equations that are sent to IMSL’s ZX3LP. 
Subroutine Fixer transforms all of the equality and 
inequality constraints into arrays that ZX3LP can use and 
calls ZX3LP. It then translates the results into horizontal 
and vertical moments at each node. 

Once this is complete the moments at each node and the 
specific loading scheme are used by BEAMCK to determine the 
maximum moment and its location on every loaded member- If 
automatic constraints was chosen 'at the beginning of 
execution then the program will automatically constrain any 
moment found to be greater than the yield moment and iterate 
until a solution is found where no moment exceeds the yield 
moment or until the maximum number of iterations is reached. 



Subroutine OUTPUT writes all grid information to a data 
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Figure 5 Macroscopic Program 
Flow Diagram 



sum 
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file named MOMENTS.DAT, All inputs are summarized and a 
ZX3LP error number is provided. Output includes all members 
ind moments, nodal loads (from decomposed general loads) and 
the maximum moment and location for each loaded member. 
Figure 5 shows the macroscopic program flow. 



SUBROUTINES 



HINT AND INTEL 

MINT initializes the nodal bending moment arrays. Any 
moment that is known to be zero because of boundary 
conditions are set to zero and all others, which are 
unknowns, are set to one. The master program GRIDS then 
calculates the number of unknowns in the problem by counting 
the non“zero elements in the array- These ones and zeros 
are used as multipliers to cause the coefficients in the 
equilibrium equations calculated in CONST to go to zero to 
comply with boundary conditions- 

INTPL sets the starting value of each nodal load to 



zero- 



ENTERL 

ENTERL is the subroutine that makes load variable 
assignments- Its operation is straight forward and requires 
little explanation- Some error traps are built i n to 
ensure all loads are in fact valid from the standpoint of 
location- There i s no trap however to ensure that the 
restrictions on the numbers of various loads are observed- 



JL 
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If these 

Arrays 
PTLOAD(i ) 
1 1 (i ) 

JJ (i ) 
LDlND(i ) 
DIVAL (i > 

DlLOC(i ) 

LD2ND(i ) 
D2AVL ( i > 

D2BVL(i > 

LD3ND ( i ) 

D3VAL ( i ) 
VMNT ( i > 
HMNT ( i > 
BHKMNT ( i ) 

BVKMNT ( i > 



restrictions are not observed the calculated 

PROGRAM VARIABLES 



The total load for node (i) 

The Y coordinate of node (i) 

The X coordinate of node (i) 

The reference node number of point load (i) 

The magnitude of point <i), a negative 
value indicates the point load is on 
vertical member, a positive value 
indicates the load is on the horizontal 
member 

The distance of point load (i ) from the 
reference node 

The reference node number of line load (i) 

The magnitude of line load (i) at the 
reference node, a positive value 
indicates the load is on the horizontal 
member 

The magnitude of line load (i) at the node 
opposite to the reference node, a positive 
value indicates the load is on the horizontal 
member 

The reference node number of distributed 
1 oad ( i ) 

The magnitude of the distributed load 

The value of the vertical moment at node (i) 

The value of the horizontal moment at node (i) 

Maximum moment on horizontal member of 
reference node <i) 

Maximum moment on vertical member of 
reference node (i) 



Table 1 Program Variables 
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BVKLOCd ) 


Location of max moment on vertical member of 
reference node <i) 


BHKLOC ( i ) 


Location of max moment on horizontal 
member of reference node (i) 


XCNST(i ) 


Location of manual constraint for 
reference node (i) 


NCNST ( i ) 


Reference node of manual constraint (i) 


FOG ( i , j ) \ 
A2 ( i , j ) \ 

Bl(i) / 

B2 ( i ) / 


These are dummy arrays for manipulating data 
for use by the linear programming subroutine 



Var i abl es 




NODESH 


Number of beam intersections in the horizontal 
direction (nodes) 


NODESV 


Number of beam intersections in the vertical 
direction (nodes) 


NODTOT 


Total number of beam intersections, all edges 
have beams along them 


VLEN 


Length of the vertical members, distance 
between nodes 


HLEN 


Length of the horizontal members, distance 
between nodes 


LPT 


Number of point loads added to the grillage 


LriLDN 


Number of linear line loads added to the 
grill age 


LDIS 


Number of distributed loads added to the 
grill age 


I TOP 


Top edge boundary condition 


ILFT 


Left edge boundary condition 



Table 1 Program Arrays and Variables 
(cont i nued ) 
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IRGT 


Right edge boundary condition 


IBTM 


Bottom edge boundary condition 


PP 


Grid load factor 


lER 


Linear programming routine error number 


lEQNM 


Number of non-zero moments in the grid 


NODEEQ 


Number of nodal equations required to 
obtain a solution 


HDRZMO 


The fully plastic moment of the horizontal 
grid members 


VERTMO 


The fully plastic moment of the vertical 
grid members 


Other 


variables used throughout the subroutines are 



dummy variables used in calculations or for keeping track of 
vector sizes and locations. 



Table 1 Program Arrays and Variables 
(continued ) 
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max i mum moment and location on each member will not be 
correct . 

There are several arrays for each load type- The total 
number of each type of load is kept track of with load 
counters- LPT is the total number of point loads entered 
and LNLDN is the number of line loads. LDIND and LD2ND 
store the reference node numbers of each point load and line 
load respectively. DIVAL and DILOC store the magnitude and 
location (distance from the reference node) of each point 
load. D2AVL and D2BVL contain the end point magnitudes of 
each line load- 

Table (2> gives a summary of load array functions. 

All array elements are indexed by load number for each 
specific type of load. Each load type has an array that 
contains the reference nodes associated with the load 
numbers. The sign of the load magnitude determines which 
member, either horizontal for positive or vertical for 
negative, the load is associated with. The reference node 
is always the node to the left or above the loaded member - 
Figure 6 shows the program flow for entering loads. 

EQN5 

Subroutine EQNS writes the nodal equations of the 
grillage and writes them to an array in a format that 
subroutine FIXER uses later to assimilate ail of the 
required data into the proper format. EQNS calculates the 
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Point Loads LDlND(i) 


Contains the reference 
node number for point load 
number i 


DIVAL <i ) * 


Contains the magnitude of 
load number i 


DlLOC(i ) 


Contains the location of 
(distance from the 
reference node) load 
number i 


Line Loads LD2ND<i) 


Contains the reference 
node number for point load 
number i 


D2AVL < i ) * 


Contains the magnitude of 
load number i at the 
reference node 


D2BVL ( i ) * 


Contains the magnitude of 
load number i at the 
opposite node 


odal Loads PTLOAD(i) 


Store the total nodal 
load (decomposed general 
loads) that the equilibrium 
equations are written from 







These values are set to the negative to indicate the load 
i s on the vertical member associated with the reference node 
and are left positive to indicate the load is on th 
hor i zontal member . 



Table 2- Load Array Function 
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Figure 6 Program Flow for Entering Loads (ENTERL) 
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coefficients of the equilibrium equations for each member. 
Since the unknowns are the end moments, a single linear 
equation that contains each unknown throughout the entire 
grillage is generated for each member. The result is a 
vector with twice the number of nodes plus one elements- 
Since each member has only two ends, all of the coefficients 
except two turn out to be zero- Only the two coefficients 
associated with the nodes that bound the member will be non- 
zero- The additional unknown is the load factor which is 
not constrained but will be maximized. The left half of 
this array represents all of the horizontal moments and 
members while the right half represents the vertical moments 
and members. Figure (7) shows the construction of the 
vector matrix for all equations written by both subroutines 
EQNS and CONST. 

In Figure <7) the first example equation represents the 
equilibrium equation to the member located between the node 
in the third and forth columns or node numbers three and 
four. The last column represents the load factor and the B 
vector represents the right hand side of the equation- 
Since the non-zero elements are located in the left half of 
the matrix we know that this equation represents a 
horizontal member. The second equation has non-zero 
elements in the right half of the matrix so it must be for a 



— “TO-- 



/ ; \ 

# 

u 

n 

k 

n 

0 
w 
n 
s 

1 

\ ! / 



c- 



A2 Matri: 



< — total # nodes — > 
1 2 3 n 

Hori zontal Section 

Example equations 

001 3 00000 
OOOOOOOOO 



-total # nodes — 



1 2 3 n 

Vertical Section 



0 0 0 0 0 0 0 0 0 
000300030 



< B > 

M 



r 

1 



Figure 7- 

vertical member. Since all of the nodes of the grillage are 
numbered from left to right we know that all of the 
horizontal members will be represented by coefficients in 
adjacent elements. Likewise, vertical members will be 
represented by elements in the right side of the matrix that 
are separated by the number of horizontal node in the 
grillage. Figure 8 shows the flow diagram for EQNS. 

CONST and C0NST2 

Subroutine CONST constructs additional constraint 
equations in order to refine and clarify the problem 
solution- Specified members are selected for these 
additional constraints. Coefficients of the equilibrium 
nodal equations for each of these specified members in the 
grillage that has unknown end moments are then calculated. 
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- 33 - 




Figure 8 EQNS Subroutine Flow Diagram 
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It has Automatic or manual options. If automatic is 
selected the program calculates the maximum moment found in 
a member as the sum of the linear combination of the 
member’s end moments as calculated from the initial nodal 
loading and the particular loading of the member. The 
coefficients of the nodal equations involving that member 
are then calculated using the location of the maximum moment 
found during this initial look at the member. In manual, the 
location of the moment is chosen by the user. In manual the 
user also chooses whether or not the extra constraint is to 
have a moment already developed (ie prior damage) or it is 
to be handled as a limit. 

C0NST2 operates exactly as subroutine EQNS to write an 
Array of Coefficients to be used by subroutine FIXER- 
Figure 9 shows the flow diagram for CONST and C0NST2- 



FIXER 

Subroutine FIXER gathers the known data and converts it 
into a form that can be used by the IMSL linear programming 
routine, ZX3LP. It writes an identity matrix that acts as 
the list of inequality constraints and it uses the arrays 
produced in EQNS and C0NST2 to write the equality 
constraints and additional constraints. All of these are 
combined in a single matrix that is used by ZX3LP. See 
Appendix B for details of this matrix. ZX3LP is then 
called. FIXER then converts the solution returned from 



- 36 - 




Figure lO Subroutine FIXER Flow Diagram 
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BEAMCP< 

Subroutine BEAMCK calculates the maximum moment and 
location for every loaded member. It accomplishes this by 
looking at each node and gathering the loads that are 
referenced to it. If no loads are referenced there will be 
no calculations for that node. It then separates the loads 
into vertical and horizontal members and assigns them to 
dummy variables. Equations for the moment as a function of 
position are available for each possible loading condition 
and position on a member. There are four possibilities. 

The linear line load may be either increasing or decreasing 
going away from the reference node and the moment must be 
calculated for either side of a point load. A constant line 
load will fir either type of linear line load equation. 

These equations are derived in Appendix C. 

The program applies the proper equation and calculates 
the moment and location. The results are stored in four 
arrays. BHKMNT and BHKLOC store the moment and location for 
the horizontal member referenced to a node and BVKMNT and 
BVKLOC store this information for the vertical members. It 
should be noted that subroutine CONST uses identical logic 
to determine which additional constraint equations to write 
when automatic constraints is selected. Figure 11 shows the 
subroutine flow for BEANCK. 

OUTPUT 



Subroutine OUTPUT consists almost entirely of write and 
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GATHER ALL THE LOADS 
ASSOCIATED WITH A NODE 
AND ASSIGN TO DUNliY VARIABLES 
IF NO LOADS SKIP THE NODE 



CHOOSE THE COR 
FOR THE LOADING CON 


RECT EQUATION 
DITION ENCOUNTERED 






USING THE END MOMENTS CALCULATED 
USING ZX3LP, CALCULATE THE MAXIMUM 
MOMENT AND LOCATION 






STORE THE MAXIMUM MOMENT AND LOCATION 
IN BHKMNT AND BHKLOC (HORZ MEMBER) 
OR BVKMNT AND BVKLOC (VERT MEMBER) 







c 



RETURN 



3 



Figure 11 BEAMCK Subroutine Flow Diagram 
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^^ormat statements that provide a hard copy of the grid 
character i sti cs, loading and the solution- A ZX3LP error 
number is also provided which can be referenced in Appendix 
B if problems are encountered- A sample output is shown in 
Figure 12- 

RUNNir4G THE PROGRAM 

Once GRIDS and its subroutines- including the IMSL 
routine ZX3LP and ZX30P (ZX30P is used by ZX3LP) , are 
installed on the operating system, the program can be run- 
It may be run interactively using keyboard input or it may 
be run using an input data file- Appendix B shows how to 
create an input data file- 

The following listing is the actual program dialogue as 
seen on the terminal when running Grids. User responses have 
been put in square brackets for cl ar i f i cat i on - The listing 
below is typical but not all inclusive- If automatic 
constraints is selected, the inputs that concern manual 
constraints is not asked for- All other inputs are similar- 

RUN GRIDSl 

PROGRAM RESTRICTED TO A MAXIMUM GRID OF 6 X 5 
IE 30 NODES TOTAL OR LESS 

THIS RESTRICTION IS IMPOSED BECAUSE OF THE MAXIMUM 
SIZE OF THE WORK VECTOR USED BY ZX3LP IS CONSTRAINED 
BY THE PROGRAMMER. 

SEE PROGRAM MANUAL FOR MORE INFO 

INPUT NUMBER OF NODES ALONG THE TOP AND SIDE OF THE GRID 
C4-41 

INPUT BOUNDARY CONDITIONS FOR THE TOP 
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LEFT. RIGHT AND BOTTOM EDGES OF THE GRID 
IN THAT ORDER 

INPUT A 1 FOR SIMPLY SUPPORTED EDGE 
INPUT A 2 FOR A CLAMPED EDGE 
ALL VALUES ARE INTEGERS 

Cl, 1, 1, 13 

INPUT HORIZONTAL MEMBER LENGTH (REAL NUM, ) 
C15. 3 

INPUT VERTICAL MEMBER LENGTH (REAL NUM.) 
C15. 3 

INPUT HORIZONTAL BEAM YIELD MOMENT (REAL) 

C 30000. 3 

INPUT VERTICAL BEAM YIELD MOMENT (REAL) 

C 30000. 3 

LOADS ARE ENTERED REFERENCED TO THE NODE 
ABOVE OR TO THE LEFT OF THE LOAD 
ALLOWED LOADING CONFIGURATIONS: 

1. POINT LOADS 

2. LINEAR LINE LOADS ALONG A MEMBER 

3. LOAD ON A PLATE ELEMENT 



ENTER THE TYPE OF LOAD TO ADD 

1 = POINT LOADS 

2 = LINEAR LINE LOAD ALONG A MEMBER 
0 = NO MORE LOADS TO ADD 

C13 

THIS ROUTINE ADDS POINT LOADS TO THE GRID. THERE MAY BE 
A MAXIMUM OF ONE POINT LOAD PER NODE PLUS ONE PER MEMBER 
THE REFERENCE NODE IS ABOVE OR TO THE LEFT OF THE LOAD 
ENTER FOUR VALUES TO DEFINE EACH LOAD 

1ST VALUE = REFERENCE NODE (INTEGER) 

2ND VALUE = VERTICAL OR HORIZONTAL MEMBER 

1 = HORIZONTAL MEMBER (INTEGER) 

2 = VERTICAL MEMBER (INTEGER) 

O = ON THE NODE (INTEGER) 

3RD VALUE = LOAD MAGNITUDE (REAL) 

4TH VALUE = DISTANCE FROM THE REFERENCE 

NODE (REAL) (IF THE 2ND VALUE = 0 
THEN THE 4TH MUST = O) 

ENTER ALL ZEROS WHEN COMPLETED WITH POINT LOADS 
C6, 1,7.5,7.53 

ENTER FOUR VALUES TO DEFINE EACH LOAD 

1ST VALUE = REFERENCE NODE (INTEGER) 

2ND VALUE = VERTICAL OR HORIZONTAL MEMBER 
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1 = HORIZONTAL MEMBER (INTEGER) 

2 = VERTICAL MEMBER (INTEGER) 

O = ON THE NODE (INTEGER) 

3RD VALUE = LOAD MAGNITUDE (REAL) 

4TH VALUE = DISTANCE FROM THE REFERENCE 
NODE (REAL) (IF THE 2ND VALUE 
THEN THE 4TH MUST = O) 

ENTER ALL ZEROS WHEN COMPLETED WITH POINT LOADS 
CO, 0,0. ,0 . 3 

ENTER THE TYPE OF LOAD TO ADD 

1 = POINT LOADS 

2 = LINEAR LINE LOAD ALONG A MEMBER 
O = NO MORE LOADS TO ADD 

CO 3 

THIS SUBROUTINE ALLOWS PLASTIC MOMENTS TO BE 
FORMED AT LOCATIONS OFF OF THE NODES OF THE 
GRILLAGE. THIS IS DONE BY WRITING ADDITIONAL 
CONSTRAINT EQUATIONS FOR THE LINEAR PROGRAMMING 
ROUTINE 

AUTO CONSTRAINTS WILL LOOK AT EACH MEMBER AND 
CONSTRAIN THE MAXIMUM MOMENT TO <= YIELD MOMENT 



MANUAL ALLOWS THE USER TO CHOOSE WHICH MEMBERS 
TO SET CONSTRAINTS FOR AND THEIR LOCATION 



INPUT A 1 FOR AUTOMATIC CONSTRAINTS 
A 2 FOR MANUAL CONSTRAINTS 
A O FOR NO CONSTRAINTS 

C2D 

CONSTRAINED MEMBERS MUST HAVE AT LEAST ONE LOAD 
ASSOCIATED WITH THEM 

INPUT THE REFERENCE NODE FOR THE ADDITIONAL 
CONSTRAINT, AND A 1 FOR HORIZONTAL MEMBER 
OR A 2 FOR VERTICAL MEMBER (INTEGERS) 

ENTER ZEROS FOR NO CONSTRAINTS 

C6. 13 

INPUT THE DISTANCE FROM THE REFERENCE NODE 
TO LOCATE THE CONSTRAINT (REAL) 

C7.53 

ENTER A 1 TO INDICATE THAT THE MEMBER HAS A 
FULLY DEVELOPED YIELD MOMENT ALREADY DEVELOPED 
NOTE: THIS TYPE MUST BE ENTERED LAST 
ENTER AO IF THE HINGE IS NOT ALREADY FORMED 
( INTEGERS) 



C03 
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At this point the program will run. All of the output is 
written to a data file that can either be written to the 
terminal or printed on the system hardcopy device. 



-43- 



GRID CHARACTERISTICS AND LOADING 

NUMBER OF NODES HORIZONTAL = 4 HORIZONTAL LENGTH = 15.00 

NUMBER OF NODES VERTICAL = 4 VERTICAL LENGTH = 15.00 

BOUNDARY CONDITIONS 



TOP 

BOTTOM 

LEFT 

RIGHT 



SIMPLE SUPPORTED 
SIMPLE SUPPORTED 
SIMPLE SUPPORTED 
SIMPLE SUPPORTED 



•HORZ BEAM PLASTIC BENDING MOMENT 



0. lOE+04 



VERT BEAM PLASTIC BENDING MOMENT 



O. lOE+04 



ZX3LP ERROR NUMBER = 0 



POINT LOADS ADDED TO THE GRID 

REF NODE DIR FM NODE DISTANCE FROM NODE MAGNITUDE LOA 

NONE USED 



LINE LOADS 



BETWEEN 


NODES 


REF NODE MAG 


OTHER NODE MAG 


LINE LOAD 


6 


AND 


7 


1.00 


2.00 


1 


6 


AND 


10 


1 . 00 


2.00 


o 

J- 


7 


AND 


11 


2.00 


1 . 00 




10 


AND 


11 


2.00 


1.00 


4 




NODE 


LOADS USED 


IN CALCULATIONS 







NODE NUMBER 
1 



NODE LOAD 
0,00 
0 . 00 
O. 00 



Figure 1 



Sample Output 
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4 


0.00 


5 


0. 00 


6 


20.00 


7 


25. 00 


8 


0.00 


9 


0.00 


10 


25.00 


11 


20.00 


12 


0.00 


13 


0. 00 


14 


0. 00 


15 


0. 00 


16 


0. 00 



ADDITIONAL CONSTRAINTS USED 

BETWEEN NODES TYPE USED 

O AND O NONE 

**** ITERATION NUMBER IS 1 ***** 
NODAL MOMENTS 



NODE NUMBER 


HORIZ MNT 


VERT MNT 




1 


0- OOE+00 


0. OOE+00 




O 

X- 


0- 00E-*"00 


0. OOE+00 






0. OOE-^00 


0. OOE+00 




4 


0. OOE+00 


0. OOE+00 




5 


0. OOE+00 


0. OOE+OO 




6 


0. lOE+04 


0. 86E+03 




7 


0. lOE+04 


0 . 1 OE+04 




8 


0. OOE+OO 


0- OOE+OO 




9 


0. OOE+00 


0- OOE+00 




10 


0. lOE+04 


0. lOE+04 




11 


0- lOE+04 


0. 86E+03 




12 


0. OOE-^00 


0. OOE+00 




13 


0. OOE+00 


0- OOE+00 




14 


0. OOE+00 


0. OOE+00 




15 


0- OOE+00 


0- OOE+00 




16 


0. OOE+00 


0. OOE+00 




THE GRID LOAD 


FACTOR = 5. 


71 




MAXIMUM MOMENTS AND 


LOCATION ON LOADED MEMBERS 




BETWEEN NODES 


MAGNITUDE 


LOCATION 




HORIZONTAL MEMBERS 


6 AND 7 


1242. 


7.905 


** 


10 AND 11 


1242. 


7.080 




VERTICAL MEMBERS 


S AND 10 


1179. 


8.985 




7 AND 11 
INDICATES > YIELD 


1179. 

MOMENT 


6.015 


** 


Figure 12 


Sample Output 


(continued ) 
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CHAPTER 3 

DISCUSSION and EXAMPLE PROBLEM 

Work done by ABBOTT (2) resulted in a program capable 
estimating grid load factors based only on nodal loads. 
These estimates turn out to be overly conser vat i ve. For any 
symmetr i cal 1 y loaded grillage one can expect the failure mode 
to also be symmetric. In applying linear programming to find 
a failure mode (or any technique that assumes nodal loading) 
a load factor and a pattern of plastic yield hinges at the 
grid intersections is developed. It should be noted that 
the solution found will not necessarily be unique. 

In maximizing the number of plastic hinges for any 
given system there is bound to be adjacent nodes with fully 
developed plastic hinges. This will be true in all cases 
except the nine node grid. If a member has two fully 
developed yield moments of the same sign at each end, then 
the moment developed within the span must be greater than at 
the ends for any consistent loading on the member- The 
implication is clear- For any consistent loading scheme 
there must be additional yield moments formed at locations 
other that at the nodes- In comparing the results of the 
program GRIDS to the results presented by Hodge (1) or 
Abbott (2) it can be seen that a program that allows for 
yield moments to be formed with in a span returns a much 
more accurate result. 

In taking a simple symmetric case and analyzing it for 
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various conditions the use of GRIDS can be demonstr at ed - 
The program's limits and a straight forward verification also 
result from this exercise. A nine node grillage was 
analyzed using a point load to compare relative load 
carrying capability and collapse modes for various boundary 
conditions and load conditions. 

The first geometry considered is the nine node grillage 
with clamped boundaries. The load was first applied at the 
node and then gradually moved toward the edge along one 
member. First, if this grillage is analyzed in the same way 
as the simple supported grid in chapter two was, the grid 
load factor is found to be 8. A yield moment forms at the 
intersection of the members in both the horizontal and 
vertical directions. With the edges clamped hinges will 
also form at the boundary end of each beam- Setting the 
work done in deflecting the node to the sum of the work done 
in forming each of the plastic hinges will give the upper 
plastic limit or collapse load of the grid. The following 
equalities are obtained. Here F is the collapse load. 

FLe = 8 Hoe (14) 



F 



SMo/L 



( 15 ) 
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Figure 13 



Collapse Mode With Plastic Hinge at Midspan 
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Now, if the point load is moved to midspan, we can 
calculate a new collapse load and allow for a plastic hinge 
to form at the point on the member that will see the maximum 
moment. Figure 14 shows the collapse mode and load 
conf iguration for a point load located half way between the 
node and the edge of the left center horizontal member. 
Because of the symmetry involved in this case it is 
relatively simple to calculate the collapse load by hand. 
Summing the total work done in deforming the grid we have 
the force applied moving through a distance of 8_L/2 and the 
center node with its effective force displaced a distance of 

F/2 e^L + F e, L/2 (16) 

And summing the work required to form the plastic moments 
gives: 



2Moe_ + 2Moe_ 6Mo8 (17) 

*.:• 2 1 

Setting (3) and (4) equal to each other and realizing that 

03 is the sum of 0. and 0^ gives: 

1 ^ 

FL(20, + 0^) = 8Mo (20 + 0^) (18) 

1 ^ 1 



F 



SMo/L 



(19) 
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The fact that this result is the same as found in (15) 
above is only coincidental and is due to the symmetry 
involved. Similar calculations can be performed for any 
position along a member. A tabular comparison of the 
results of the program Grids and the results of a program 
that allows only nodal loading as a function of location of 
the point load is presented in Table 3- For the application 
of a point load and no intermediate yield moment, it is 
readily apparent that the farther the load moves from the 
node the more error is introduced in computing the collapse 
load. This occurs because the yield moment that must form 
at the point load will be formed with much less load as the 
load moves away from the node. 



LOCATION * 
IN 7. OF 
LENGTH 



CLAMPED NINE NODE LOAD FACTORS 

LOAD FACTOR WITH LOAD FACTOR WITHOUT 

PLASTIC MOMENTS PLASTIC MOMENTS 

IN THE MEMBER IN THE MEMBER 



10 

30 

50 

70 

90 

* measured from 



767.68 

747.25 

800. 

952.38 
2222 . 

the center node 



888.89 
1142.86 
1600. 
2666. 67 
8000. 



Table 3 Load Factor for Various Load Locations 
For Plastic Moments in a Member and For 
no Plastic Moments Except at nodes 
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If the load is normalized to a total nodal load ot one 
for each position along the member, the relative load 
carrying capability of the grillage can be plotted as a 
ratio of the load factor at a point divided by the load 
factor for nodal loads only versus position along a member. 
Table 4 lists the data for the nine node example and Figure 
14 shows the result graphically. 



YIELD HINGE IN MEMBER 



PP 

PP 



n 



LOCATION (7.) 
FROM CENTER 
NODE 



NO YIELD HINGE IN MEMBER 
PP LOCATION (7.) 

PP FROM CENTER 

^ NODE 



1 


0 


0.8636 


10 


0.6538 


30 


0.5 


50 


0.3571 


70 


0.2778 


90 



1 

1 

1 

1 

1 

1 



O 

10 

30 

50 

70 

90 



TABLE 4 Normalized Load Factor VS Position of Hinge 
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Figure 14 Normalized Collapse Load VS Position of Hinge 
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The second type of problem tested was a grillage that is 
a more realistic model of a ship's side- A 6 X 5 grillage 
with a line load was analyzed. A line load is a reasonable 
approximation of a ridge of ice along a ship's side- The 
line load was placed along a series of horizontal members 
and the load factor and failure mode were calculated- The 
load was moved vertically downward and the calculations were 
repeated for different positions along the vertical beam- 
Note that once the line load is off of the horizontal 
member, it acts like a series of point loads along a line 
across each of the vertical members- Figure 15 shows the 
loading and a typical failure mode encountered- Table 5 
shows the comparison of the load factors given by grids to 
those that are the result of a program that allows no 
intermediate hinges to form- It is interesting to note that 
the load factor remains more nearly constant when 
intermediate yield moments are allowed to form- 

The final problem to be considered is one that was 
considered by Abbott (2). Here we must convert the 
distributive load into line loads on the members- The 
problem is an example that illustrates how GRIDS might be 
used in considering the strength of a ship's side that is 
anticipated to operate in a particular ice environment- For 
this problem consider that the ship is to operate in a 
uniform ice sheet five feet thick- Further assume that the 
ice has a local crushing strength of 300 psi and that 



LOCATION (7.) 



WITH INTERMEDIATE 
PLASTIC MOMENTS 



NO INTERMEDIATE 
PLASTIC MOMENTS 



o 


384.00 


400. 00 


10 


380 . w6 


413. 79 


30 


353.28 


444.44 


50 


337. 78 


444.44 


70 


330. 78 


392. 16 


90 


330.83 


350.88 



Table 5 Load Factor for 5X6 VS Position of Line Load 
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Figure 15 Failure Mode for Line Load on Horizontal Member 
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HORIZONTAL MEMBERS 




4 




VERTICAL MEMBERS 



Failure Mode and Deformation Profile for 
with a Horizontal Line Load Between 
Horizontal Members With Clamped Boundaries 



Figure 16 



6X5 
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crushing is the predominant failure mode. Analyze a ship 
which has the following character i st i cs- 



Frame Spacing 
Bulkhead Spacing 
Distance Between Decks 
Number of longitudinals 
Stringer Modulus 
Frame Modulus 
Yield strength 



21 inches 
7 feet 
15 feet 

2 (evenly spaced) 
.5 in" 

-T 

.8 in“ 

50000 psi 



Figure 17 shows the conf i gurati on for this problem- 



KE 21.^4^ 21 - rX: 21 21 




CZ] 



ICE LOADED 



Figure 17 Ice Loaded Ship's Side 
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Converting the distributive ice loads on this grillage 
we get lineloads around each plate element- Each plate has a 
total pressure of 3.78 X lO'^ pounds on it. Assuming this 
reduces to a uniform line load around the perimeter of each 
plate results in a line load of 583.33 pounds per inch. 
Entering the data and selecting automatic constraints 
results gives a Load Factor of 0.175. With this load factor 
we know that this grillage can support less than twenty 
percent of the anticipated load. The structure must be 
strengthened. Strengthening can be accomplished in several 
ways- Increasing plate thickness, frame of longitudinal 
member sizes or using stronger materials are all viable 
options that will result in an increased maximum yield 
moment of each member. Increasing the yield moment of each 
member results in a grillage that will support the applied 
1 oad. 

Additional insight can ge gained by varying the strength 
of either the horizontal or vertical members alone and 
observe the structural behavior. With the same load applied 
as in the problem above, the horizontal yield moment was set 
at a value of 100,000 and the vertical bending moment was 
varied through a range of values of 20,000 to 90,000. All 
of these values give the same Load Factor of O. IS- This 
suggests that for the applied loading the grillage is much 
weaker in the horizontal direction, which makes sense when 
the geometry is considered. In looking at the failure mode 
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provided by GRIDS we observe that plastic hinges are 
indicated for all boundary nodes- We also know that there 
are a number of yield moments given that do not actually 
occur, due to the degree of indeterminancy and linear 
programming methods (Chapter 4 contains additional 
discussion on this point). We can conclude from the above 
that the moments in the vertical members are the ones that 
are not actually formed. 

Now, if we hold the yield moment of the vertical 
member constant and vary the yield moment of the horizontal 
member we can observe the grillage behavior as a function of 
the horizontal member size. Figure 16 shows how the load 
factor varies as the horizontal member yield moment for this 
example. 




Load Factor 
4 

Vertical Yield Moment = 90,000 in 
Figure IS Load Factor Versus Horizontal Member Strength 
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CHAPTER 4 

SUMMARY OF RESULTS 

As illustrated in Chapter 3 the results that Grids 
gives are significantly more accurate than for a program that 
allows only nodal loading. In choosing to use linear 
programming to solve a problem with intermediate yield 
hinges, an iterative approach is also required. If one 
were to write constraint equations for forming the plastic 
hinges within a member directly, as a function of a general 
combination of point and line loads, the resulting equations 
would be non-linear. In order to use linear programming and 
still allow plastic hinges in a member, the end moments must 
be used. So, the problem is first solved using only nodal 
loads to calculate the end moments. Constraint equations 
are then written that are a linear combination of loading 
and end moments. The assumption is that the final end 
moments will not change the basic shape of the moment 
distribution. The end moment only changes the slope of the 
line that the moment distribution lies on. From this one 
can see that if the end moments, as given by GRIDS, for a 
particular member are equal then the location of the maximum 
moment is exact. Hence the position of the yield moment in 
that member is correct. 

For point loads only the position of the failure is 
known and the program solution is accurate. In the case 
of line loads GRIDS gives a "linearized" approx i mat i on . 
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Calculating the approximate error can be accomplished by 
comparing the location of the maximum moment based on the 
moment distribution and member end moments to the position 
that the additional constraint was written for. Doing this 
assumes that the yield moment that forms at the location of 
the constraint does not change the moment distribution. For 
a constant line load, the location error in percent of the 
length of the member is approximately one half of the 
percent difference between end moments. In other words, if 
the end moments are within twenty percent of each other then 
the plastic yield moment will be within ten percent of the 
location given by GRIDS. Further more, it will be shifted 
toward the end of the member with the larger moment - 

The location that GRIDS calculates for the maximum will 
be exact for point loads and constant line loads on a 
sysmmetrical grid with symmetrical loading. In non- 
symmetrical grids with multiple load types the locations of 
the maximum moments are approximate but very close to 
actual- This occurs because of the way in which linear 
programming calculates the end moments of each member. If 
the end moments of a member are unchanged between the first 
and second iterations as labeled in the output, then the 
location of the maximum moment in that member calculated by 
grids is exact- If one of the end moments of the member is 
reduced then the effect is to lower that end of the moment 



curve until the maximum moment as constrained is found - 



The 
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curve keeps its original shape however. The error 
introduced is small as long as the end point magnitudes of a 
member's linear line loads do not differ more than fifty 
percent . 

GRIDS offers both automatic and manual modes for 
developing the additional constraint equations that actually 
add the plastic yield moments in a member. This gives added 
flexibility in evaluating a particular grillage. By 
choosing the manual mode and choosing equality constraints a 
failure can be force to form at a particular location. By 
forcing a failure, or by having an already formed plastic 
hinge in a member, the load carrying capability of a 
partially failed grid can be examined. This technique may 
be desirable to evaluate a grillage that has undergone a 
local impact on a member so that a hinge may be formed at a 
location first that might otherwise not reach its yield 
moment until much higher general loads are applied. If, 
however, a plastic hinge is placed in a member with no loads 
associated with it the output is meaningless and a ZX3LP 
linear programming error will occur - 

As discussed previously, the failure mode provided by 
GRIDS is not unique. All of the systems encountered contain 
more unknowns than equations and therefore cannot be solved 
directly. The way in which this linear system is solved and 
maximized causes selected inequality constraints to be 
evaluated as equality constraints. The number used in this 



- 60 - 



fashion is equal to the degree of i ndetermi nancy of the 
system. This means that not all of the locations that are 
listed in the output may actually reach the yield moment. 
But, the output provided is a maximized and feasible 
solution. In other words, there may be fewer yield moments 
formed in an actual failure but the total work done is 
correct. This could be accounted for by different 
magnitudes of rotation of the hinges that do form- 
FUTURE WORK 

GIRDS can be modified to handle a much broader base of 
problem than those which have been considered up to now. 
There are some modifications and additions that could be 
made to enhance the program's performance. 

The most useful feature that could be added would be 
the capability to enter distributed loads. This task should 
not be too simplistic, however. Distributed plate loads do 
not transfer to the edges of the plate as linear line loads 
except in the case of circular plates and should not be 
handled in this way. One might approach this problem by 
linearizing or approxi mating the Fourier series solution for 
loaded plates. 

Another area that improvement could be made is to allow 
for different yield moments in particular members. Adding 
arrays to keep track of this specific information and 
calling new variables as appropriate would be tedious but 
not particularly difficult from a programming standpoint. 
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APPENDIX A 



PROGRAM LISTING AND DIRECTIONS 
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Grids is a very simple program to use- In the 
interactive mode it is self explanitory and requires no 
further information to run. To run the program more quickly 
when repetitive inputs are desired, an input data file may be 
used- Figure 19 shows the construction of such a file. One 
need only assign that file to be the source of input to run 
the program with all input read from the input file without 
further action from the user. 

The Grillage size limits placed on GRIDS was 
primarily for the convinence of the programmer. Arrays may 
be redi mensi oned to solve larger grids. The criteria 
provided in Appendix B for the IMSL linear programming 
routine must be met , however - 
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SAMPLE INPUT DATA FILE (WITH EXPANATIONS) 



5,4 


GRID DIMENSIONS ON THIS LINE 


o O O 

« wC. 9 X. 4 wC. 


BOUNDARY CONDITIONS 


21. 


HORIONTAL LENGTH 


60. 


VERTICAL LENGTH 


30000. 


HORIZONTAL YIELD MOMENT 


90000 . 


VERTICAL YIELD MOMENT 


o 


TYPE OF LOAD (1,2,0) 


6, 1,53.3,53.3\ 





6,2,53.3,53.-3 \ 



7.2.53.3.53.3 
7 1 ^^3 3 ‘=i3 3 
8, 1,53.3,53.3 

8.2.53.3.53.3 

Q ^ *r -r 

r jC- n ■ *_* 4 • *_* 

Q 1 «=:-r “r 't 

10.2.53.3.53.3 

11. 1.53.3.53.3 

12. 1.53.3.53.3 

13. 1.53.3.53.3 
14, 1,53.3,53.3/ 


\ 

\ 

REFERENCE NODE, DIRECTION (1.2), 
VALUE FOR POINT LOAD, OR A-NODE 
VALUE FOR LINE LOAD, LOCATION FOR 
/ POINT LOAD, OR B-NODE VALUE FOR LINE 

/ 


0,0,0. ,0. 


NO MORE LOADS OF THIS TYPE TO ADD 


0 


NO MORE LOADS TO ADD 


1 


TYPE OF CONSTRAINTS (1,2,0) 



There may be no blank lines in the data -file. Lines were left 
blank for cl ar i f i cat i on . There are four additional lines required 
if manual constraints are used. 



Figure 19 Sample Input Data File 



no n on 
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C THIS BLOCK IS THE COMMON MEMORY SHARED BY ALL SUBROUTINES 
C IT IS INCORPORATED BY USING THE INCLUDE ATATEMENT 

COMMON /A/ PTLOAD (30) , NODTOT, NODESH, NODESV, 

+ 11(30) , JJ (30) , VLEN,HLEN, XDCNST(30) , 

+ LD1ND(60) ,D1VAL(60) ,D1L0C(60) , LPT, LD2ND (60) , 

+ D2AVL(60) ,D2BVL(60) , LNLDN, LDIS, NDCNST (30) , 

+ LD3ND (60) , D3VAL (60) , I AUTO, NEQST 

COMMON /B/ VMNT(30) ,HMNT(30) jHORZMO, VERTMO 
COMMON /C/ ITOP, ILFT, IRGT, IBTM, ITERAT 
COMMON /D/ PP, lER, lEQNM, NODEEQ,NITRAT 
COMMON /E/ F0G(12,61) ,A2(31,61) ,B1 (12) ,B2(31) 

COMMON /F/ LIN(40) ,MN, IA,LMAX,M1 ,M2, IWl, IRW1,LP2 

COMMON /GGG/ BHKMNT (30) , BHKLOC (30) , BVKMNT (30) , BVKLOC (30) 

C THIS PROGRAM SOLVES THE PROBLEM OF FINDING THE MAXIMUM LOAD FACTOR 

C TO BE APPLIED TO A PARTICULAR LOADING CONDITION OF A RECTANGULAR Gf 

C IT SOLVES THE PROBLEM BY CONVERTING THE RESULTING SET OF EQUATIONS 

C INTO A LINEAR PROGRAMMING PROBLEM. AT THE PRESENT TIME THE MAXIMUr 

C SIZE OF GRID THAT CAN BE HANDLED IS A 6 X 5 GRID (IE 30 NODES ) 

C THIS CONTRA I NT IS ARTIFICIAL AND IS DUE SOLELY TO THE VERY LARGE 

C WORK VECTOR REQUIRED BY THE LP SUBROUTINE ZX3LP 

C SEE IMSL LIBRARY REFERENCE MANUAL #8 FOR FURTHER INFO. 

DIMENSION A (83, 39) ,B (83) ,C(39) ,PS0L(83) ,DS0L(83) , 

+ RW (7200) , IW (240) 



FOR' 



0PEN(UNIT=15,FILE='M0MENTS.DAT' , STATUS=' NEW' ) 

1 PRINT*, 'PROGRAM RESTRICTED TO A MAXIMUM GRID OF 6 X 5' 

PRINT*, 'IE 30 NODES TOTAL OR LESS' 

PRINT*, 'THIS RESTRICTION IS IMPOSED BECAUSE OF THE MAXIMUM' 
PRINT*, 'SIZE OF THE WORK VECTOR USED BY ZX3LP IS CONSTRAINED' 
PRINT*, 'BY THE PROGRAMMMER. ' 

PRINT*, 'SEE PROGRAM MANUAL FOR MORE INFO' 

INPUT SIZE OF GRID TO BE SOLVED 
PRINT*,' 

PRINT*, ' 

PRINT*, ' INPUT NUMBER OF NODES ALONG THE TOP AND SIDE OF THE GRID' 
READ (5,*) NODESH, NODESV 

NODTOT=NODESH*NODESV 
IF (NODTOT. GT. 30) GOTO 1 

INPUTCONDITIONS 



noon naan hjn 
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PRINT*, ' INPUT BOUNDARY CONDITIONS FOR THE TOP' 
PRINT*, 'LEFT, RIGHT AND BOTTOM EDGES OF THE GRID' 
PRINT*,' IN THAT ORDER' 

PRINT*, ' INPUT A 1 FOR SIMPLY SUPORTED EDGE' 

PRINT*, ' INPUT A 2 FOR A CLAMPED EDGE' 

PRINT*, 'ALL VALUES ARE INTEGERS' 

PRINT*, ' 

PRINT*,' 

READ(5,*) ITOP, ILFT, IRGT, IBTM 

I=ITOP 

J=ILFT 

K=IRGT 

L=IBTM 

IFd.EQ. l.OR. I.EQ.2. AND. J.EQ. l.OR. J.EQ.2. 

+AND.K.EQ. 1.0R.K.EQ.2. AND.L.EQ. l.OR.L. 

PRINT*, 'BOUNDARY CONDITIONS NOT ENTERED CORRECTLY.' 
PRINT*,' TRY AGAIN' 

GO TO 2 
CONTINUE 

PRINT*, ' INPUT HORIZONTAL MEMBER LENGTH (REAL NUM. > ' 
PRINT*, ' 

READ (5,*) HLEN 

PRINT*, ' INPUT VERTICAL MEMBER LENGHT (REAL NUM.)' 
PRINT*, ' 

READ (5,*) VLEN 

PRINT*, ' INPUT HORIZONTAL BEAM YIELD MOMENT (REAL)' 
PRINT*, ' 

READ(5,*> HORZMO 

PRINT*, ' INPUT VERTICAL BEAM YIELD MOMENT (REAL) ' 
PRINT*, ' 

READ (5,*) VERTMO 
SET LOAD COUNTERS TO ZERO 
LPT=0 
LNLDN=0 
LDIS=0 
IEQNM=0 

NOW TO INITIALILZE THE LOAD AND MOMENTS ARRAYS 
CALL MINT 



CALL INPTL 

NOW TO ENTER ALL THE LOADS APPLIED TO THE GRID 



CALL ENTERL 



NOW TO ENTER ADDITIONAL CONSTRAINT INEQUALITIES 



n n n cn on 
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SEE USER'S HANDBOOK FOR FURTHER INFO 

PRINT*, 'THIS SUBROUTINE ALLOWS PLASTIC MOMENTS TO BE' 
PRINT*, ' FORMED AT LOCATIONS OFF OF THE NODES OF THE' 
PRINT*, 'GRILLAGE. THIS IS DONE BY WRITING ADDITIONAL' 
PRINT*, 'CONSTRAINT EQUATIONS FOR THE LINEAR PROGRAMMNIG' 
PRINT*, 'ROUTINE' 

PRINT*, 'AUTO CONSTRAINTS WILL LOOK AT EACH MEMBER AND' 
PRINT*, ' CONSTRAIN THE MAXIMUM MOMENT TO <= YIELD MOMENT' 
PRINT*, ' 

PRINT*, ' 

PRINT*,' MANUAL ALLOWS THE USER TO CHOOSE WHICH MEMBERS' 
PRINT*,' TO SET CONSTRAINTS FOR AND THEIR LOCATION ’ 
PRINT*,' 

PRINT*,' 

PRINT*, ' INPUT A 1 FOR AUTOMATIC CONSTRAINTS' 

PRINT*,' A 2 FOR MANUAL CONSTRAINTS' 

PRINT*,' A O FOR NO CONSTRAINTS 

READ (5,*) lAUTO 
IF(IAUT0.EQ.2) THEN 
CALL CONST 
END IF 



CONTINUE 

NOW TO CALCULATE THE NUMBER OF NODAL EQUATIONS 

NODEEQ= (NODESH-2) * (NODESV-2) 

C 

C NOW TO CALCULATE THE NODAL EQUATION OF THE GRID AND WRITE THEM TO 

C THE MATRIX FOG 

C 

CALL EQNS 



C WE MUST CALCUALTE THE DIMENSIONS OF ALL ARRAYS BEFORE 

C FIXER BECAUSE ZX3LP REQUIRES IT 

C NOW TO CALCULATE AND TRACK THE TOTAL NUMBER OF VARIABLES 

C IE TOTAL NUMBER OF NO-ZERO MOMENTS + LOAD FACTOR 

C 

L=1 

DO 15, I=l,NODTOT 

IF(HMNT(I) .LT.0,5) GO TO 15 

LIN(L)=I 

L=L*1 

15 CONTINUE 

LP2=L-1 

DO 16 I=l,NODTOT 

IF<VMNT (I) .LT.0.5) GO TO 16 
LIN(L)=I 



noon 
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L=L+1 

16 CONTINUE 

LMAX=L 

NOW TO MASSAGE THE INFORMATION CALCULATED AND STORED BY 
THE PROGRAM SO THAT IT CAN BE USED BY ZX3LP 

M1=LMAX+IEQNM-1 
M2=N0DEEQ 
IA=Ml+M2+2 

I RW 1 = I A* I A+3*M 1 +2*M2+4 
IWl=2*M2+3*Ml+4 
MN=M1+M2 

C ADJUST Ml AND M2 IF THERE ARE EQUALITY CONSTRAINTS 
M1=M1-NEQST 
M2=M2+NEQST 

CALL FIXER ( A, B, C, PSOL, DSOL, RW, IW) 

C CALL BEAM CHECKING ROUTINE 
CALL BEAMCK 

ITERAT=ITERAT+1 
C CALL OUTPUT 

CALL OUTPUT 

BMAX=0.0 

DO 999 I=l,NODTOT 

I F ( BHKMNT ( I ) . GT . HOR ZMO ) THEN 
BMAX=BHKMNT(I) 

END IF 

I F ( BVKMNT ( I ) . GT . VERTMO ) THEN 
BMAX=BVKMNT ( I ) 

ENDIF 

999 CONTINUE 

C LOOP TO WRITE EXTRA CONSTRAINT EQUATIONS IF AUTOMATIC IS SELECTED 
IF ( lAUTO. EQ. 1 . AND. BMAX. GT. O. O. AND. 

+ ITERAT.lt. 2) THEN 

CALL CONST 
GO TO 5 
ENDIF 

END 



SUBROUTINE ENTERL 
INCLUDE ’COMMON. FOR' 



PRINT ’ LOADS ARE ENTERED REFERENCED TO THE NODE 
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PRINT *, ' 
PRINT ♦, ' 
PRINT *, ' 
PRINT ♦, ’ 
PRINT *, ' 
PRINT *, ’ 



ABOVE OR TO THE LEFT OF THE LOAD’ 
ALLOWED LOADING CONFIGURATIONS:’ 

1. POINT LOADS’ 

2. LINEAR LINE LOADS ALONG A MEMBER’ 

3. LOAD ON A PLATE ELEMENT’ 



90 



PRINT *, ’ 
PRINT *, ’ 
PRINT *, ’ 
PRINT *, ’ 
PRINT ’ 
PRINT *, ’ 



ENTER THE TYPE OF LOAD TO ADD’ 

1 = POINT LOADS’ 

2 = LINEAR LINE LOAD ALONG A MEMBER’ 
O = NO MORE LOADS TO ADD’ 



READ (5,*) MM 



IF 


(MM ,E0. 


0) 


GO 


TO 


150 


IF 


(MM .EQ. 


1) 


GO 


TO 


100 


IF 


(MM .EQ. 


2) 


60 


TO 


110 



100 CALL PNTLD (PTLOAD, NODTOT, NODESH, NODESV, VLEN 
+ ,HLEN, II, JJ,LD1ND,D1VAL,D1L0C,LPT) 

GO TO 90 

110 CALL LNLD (PTLOAD, NODTOT, NODESH, NODESV, VLEN 
+ ,HLEN, II, JJ,LD2ND,D2AVL,D2BVL,LNLDN) 

60 TO 90 



150 RETURN 
END 



SUBROUTINE PNTLD (PTLOAD, NODTOT, NODESH, NODESV, VLEN 
+ ,HLEN, II, JJ,LD1ND,D1VAL,D1L0C,LPT) 

DIMENSION PTLOAD (NODTOT) ,D1L0C(N0DT0T*2) 

DIMENSION D1VAL(NODTOT*2) 

DIMENSION 1 1 (NODTOT) , JJ (NODTOT) , LDIND (N0DT0T*2) 



PRINT *,’THIS ROUTINE ADDS POINT LOADS TO THE GRID. 

+ THERE MAY BE’ 

PRINT *,’A MAXIMUM OF ONE POINT LOAD PER NODE PLUS ONE 
+PER MEMBER’ 

PRINT *,’THE REFERENCE NODE IS ABOVE OR TO THE LEFT OF 
+ THE LOAD’ 

PRINT *, ’ENTER FOUR VALUES TO DEFINE EACH LOAD’ 

PRINT *,’ 1ST VALUE = REFERENCE NODE (INTEGER)’ 

PRINT ’ 2ND VALUE = VERTICAL OR HORIZONTAL MEMBER 



f 




1 
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PRINT 


* ' 

1 


PRINT 


* ’ 


PRINT 


* ' 


PRINT 




PRINT 




PRINT 


* ' 

* 


PRINT 





1 

O 

3RD VALUE = 
4TH VALUE = 



= HORIZONTAL MEMBER (INTEGER)’ 
= VERTICAL MEMBER (INTEGER)’ 

= ON THE NODE (INTEGER)’ 

LOAD MAGNITUDE (REAL) ’ 

DISTANCE FROM THE REFERENCE’ 
NODE (REAL) (IF THE 2ND VALUE = 
THEN THE 4TH MUST = 0) ’ 



0 ’ 



PRINT ’ENTER ALL ZEROS WHEN COMPLETED WITH POINT LOADS’ 



READ (5,*) NDNUM, IDIR, VAL.RLOCA 



IF (NDNUM .EQ. O) GO TO 130 
IF (NDNUM .GT. NODTOT) GO TO 100 

IF (II (NDNUM) .EQ. NODESV AND IDIR .EQ. 2) GOTO 110 
IF (JJ (NDNUM) .EQ. NODESH AND IDIR .EQ. 1) GOTO 110 
IF (IDIR .EQ. 2 .AND.RLOCA .GT. HLEN) GO TO 120 
IF (IDIR .EQ. 1 .AND. RLOCA .GT. VLEN) GO TO 120 
IF (IDIR .EQ. 2) THEN 

RLOCA = -1.0 * RLOCA 
END IF 

LPT = LPT + 1 

C NEGATIVE RLOCA I NL I CATES LOAD ON THE VERTICAL MEMBER 

LDIND(LPT) = NDNUM 
DIVAL(LPT) = VAL 
DILOC(LPT) = RLOCA 

C ADD LOAD IF IT IS ON THE NODE 

IF (IDIR .EQ. O ) THEN 

PTLOAD ( NDNUM ) = PTLOAD ( NDNUM ) +VAL 
GO TO 5 
END IF 

C SET DUMMY VARIABLES FOR VERT OR HORIZ MEMBER AS APPROP 

IF (IDIR .EQ. 1) THEN 
NB=NDNUM+1 
RLEN = HLEN 
ELSE 

NB=NDNUM+NODESH 
RLEN = VLEN 
END IF 

C CALCULATE AND ADD NODE LOADS FOR POINT LOAD APPLIED 

PTLOAD (NB) = PTLOAD (NB) + (ABS (RLOCA ) *VAL/RLEN) 

PTLOAD (NDNUM) = PTLOAD (NDNUM) + (RLEN-ABS (RLOCA) ) *VAL/RLEN 

GO TO 5 

100 PRINT *,'THE LAST NODE ENTERED IS NOT ON THE GRID,’ 
PRINT *, ’RE-ENTER ALL DATA FOR THE LAST LOAD’ 

GO TO 5 
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no PRINT *,'THE LAST LOAD ENTERED IS NOT ON THE GRID, ^ 
PRINT RE-ENTER ALL DATA FOR THE LAST LOAD' 

GO TO 5 



120 PRINT *, 'LOAD LOCATION IS NOT ON THE MEMBER ADJACENT TO' 
PRINT *,'THE REFERENCE NODE, RE-ENTER ALL DATA FOR THE' 
PRINT *,'LAST LOAD' 

GO TO 5 



130 RETURN 
END 



SUBROUT I NE LNLD ( PTLOAD , NODTOT , NODESH , NODES V , VLEN 
+ ,HLEN, II, JJ,LD2ND,D2AVL,D2BVL,LNLDN) 

DIMENSION PTLOAD (NODTOT) , D2AVL (NODTOT) , D2BVL (NODTOT) 

DIMENSION LD2ND (NODTOT) , II (NODTOT) , JJ (NODTOT) 



PRINT *, ’ THIS ROUTINE ADDS LINEAR LOADS, THE LOADS MUST BE 
+0N THE' 

PRINT GRID MEMBER. THE MAXIMUM NUMBER OF LINEAR LOADS 

+ ALLOWED' 

PRINT IS ONE PER GRID MEMBER. THE REFERENCE NODE IS 

+ ABOVE OR' 

PRINT ' TO THE LEFT OF THE LOAD. 



5 



PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 



*, 'ENTER FOUR VALUES TO DEFINE EACH LOAD 

1ST VALUE = REFERENCE NODE (INTEGER) 

2ND VALUE = VERTICAL OR HORIZONTAL MEMBER 
' 1 = HORIZONTAL (INTEGER) 

' 2 = VERTICAL (INTEGER) 

3RD VALUE = MAGNITUDE AT REFERENCE NODE (REAL)' 
' 4TH VALUE = MAGNITUDE AT OPPOSITE NODE (REAL) ' 



PRINT ♦, 'ENTER ZEROS FOR ALL FOUR VALUES WHEN COMPLETED 
READ (5,*) NDNUM, IDIR, AVAL,BVAL 



IF (NDNUM .EQ. 0) GO TO 130 
IF (NDNUM .GT. NODTOT) GO TO 100 

IF (II (NDNUM) .EQ. NODESV AND IDIR .EQ. 2) GOTO 110 
IF (JJ (NDNUM) .EQ. NODESH AND IDIR .EQ. 1) GOTO 110 
IF (IDIR .EQ. 2) THEN 
BVAL = -1.0*BVAL 
AVAL=-1 . 0*AVAL 
END IF 

LNLDN = LNLDN + 1 

NEGATIVE BVAL OR AVAL INDICATES LOAD ON THE VERTICAL MEMBER 



C 
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LD2ND(LNLDN) = NDNUM 
D2AVL(LNLDN) = AVAL 
D2BVL(LNLDN) = BVAL 



C SET DUMMY VARIABLES FOR VERT OR HORIZ MEMBER AS APPROP 

IF (IDIR .ED. 1) THEN 
NB=NDNUM+1 
RLEN = HLEN 
ELSE 

NB=NDNUM+NODESH 
RLEN = VLEN 
END IF 



BVAL=ABS(BVAL) 

AVAL=ABS<AVAL) 



IF (AVAL .GT. 


BVAL) THEN 


AA 


= 6.0 




BB 


= 3.0 




W = 


BVAL 




P = 


AVAL - 


BVAL 


ELSE 






AA 


= 3.0 




BB 


= 6.0 




W = 


AVAL 




P = 


BVAL - 


AVAL 



ENDIF 



C CALCULATE AND ADD NODE LOADS FOR LINE LOAD APPLIED 

PTLOAD (NDNUM) = PTLOAD (NDNUM) + ( ( W»RLEN/2) +P*RLEN/BB) 
PTLOAD(NB) = PTLOAD (NB) + ( (W»RLEN/2) + P*RLEN/AA) 



GO TO 5 



100 



PRINT *,'THE LAST 
PRINT 'RE-ENTER 
GO TO 5 



NODE ENTERED IS NOT ON THE GRID,' 
ALL DATA FOR THE LAST LOAD' 



110 



PRINT *,'THE LAST 
PRINT 'RE-ENTER 
GO TO 5 



LOAD ENTERED IS NOT ON THE GRID, ’ 
ALL DATA FOR THE LAST LOAD' 



130 RETURN 
END 



SUBROUTINE MINT 



INCLUDE 'COMMON. FOR' 



n n n n 
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150 

100 

200 

300 

10 

11 



SUBROUTINE TO INITIALIZE THE BENDING MOMENT ARRAYS BASED ON THE 
GIVEN BOUNDARY CONDITIONS 

K=0 

DO 11 I=1,N0DESV 

DO 10 J=1,N0DESH 
NDNUM=K*NODESH + J 

IFdTOP.EQ. l.AND. I.EQ. 1) GO TO lOO 
IF(IT0P.EQ.2. AND. I.ED. 1) GO TO 150 
IFdBTM.EQ. l.AND. I.EQ. NODESV) GO TO 100 
IFdBTM.EQ.2. AND. I.EQ. NODESV) GO TO 150 
IF (ILFT.EQ. l.AND. J.EQ. 1) GO TO 100 
IFdLFT.EQ.2. AND. J.EQ. 1) GO TO 300 
IF (IRGT.EQ. l.AND. J.EQ. NODESH) 60 TO 100 
IFdRGT.EO.2. AND. J.EQ. NODESH) GO TO 300 
VMNT<NDNUM)=1.0 
HMNT(NDNUM)=1.0 
II (NDNUM)=I 
JJ (NDNUM)=J 
GO TO 10 

IF ( J. EQ. 1. OR. J.EQ. NODESH) THEN 
60 TO 100 

ELSE 

GO TO 200 

END IF 

VMNT<NDNUM)=0.0 
HMNT(NDNUM)=O.0 
II (NDNUM)=I 
JJ (NDNUM)=J 
GO TO 10 

VMNT(NDNUM)=1.0 
HMNT<NDNUM)=0.0 
II (NDNUM)=I 
JJ (NDNUM)=J 
60 TO lO 

VMNT(NDNUM)=0.0 
HMNT(NDNUM)=1.0 
II (NDNUM)=I 
JJ (NDNUM)=J 
CONTINUE 

K=K+1 

CONTINUE 

RETURN 

END 

SUBROUTINE INPTL 



INCLUDE 'COMMON. FOR 



non M non 



-74- 



C THIS ROUTINE INITIALIZES ALL NODE LOADS TO ZERO 

DO 5 J=l,NODTOT 
PTL0AD(J)=0.0 
5 CONTINUE 
RETURN 
END 

C 

C THE FOLLOWING SUBROUTINE MASSA6ES THE INFORMATION CALCUALTED AND 

C STORED BY THE MAIN PROGRAM UNITS SO THAT IT CAN BE PASSED TO THE 

C LINEAR PROGRAMMING SUBROUTINE ZX3LP 

C 
C 

SUBROUTINE FIXER ( A, B, C, PSOL, DSOL, RW, IW) 

DIMENSION AdA.LMAX) ,B<IA) ,C(LMAX) ,PSOL(MN) , DSOL < I A) , 

+ RW(IRWl) , IW(IWl) 



0 



INCLUDE ’COMMON. FOR’ 

PART I — WRITE THE NODAL MOMENTS INEQUALITY CONSTRAINT EQNS 

DO 10 I=1,LMAX-1 

DO 20 J=1,LMAX 

IF(I.EQ.J) THEN 

A(I, J)=l. 

ELSE 

A(I , J)=0. 

END IF 

CONTINUE 

IF(I.LE.LP2) THEN 

B(I)=H0RZM0*2. 

ELSE 

B(I)=VERTM0*2. 

ENDIF 

CONTINUE 



NEXT WRITE THE ADDITIONAL CONSTRAINT EQUATIONS AS SELECTED BY THE 
USER IN THE LOAD ENTERING SEQUENCE 
IF(IEQNM.EQ.O) GO TO lOO 
K=LMAX+IEQNM-1 
DO 30 J=1,LMAX 
LL=1 

IF(J.EQ.LMAX) THEN 
L=2*NQDT0T+1 
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ELSE IF(J.LE.LP2) THEN 
L=LIN(J) 

ELSE 

L=LIN( J) +NODTOT 
END IF 

DO 40 I=LMAX,K 

A(I, J)=A2(LL,L) 
LL=LL+1 

40 CONTINUE 

30 CONTINUE 



LL=1 

DO 50 I=LMAX,K 

B(I)=B2(LL) 

LL=LL+1 

50 CONTINUE 

100 CONTINUE 

C 

C NOW TO WRITE THE NODAL EQUATIONS TO MATRIX A 

C 

K=LMAX+IEQNM 
DO 60 J=1,LMAX 
LL=1 

IF(J.EQ.LMAX) THEN 
L=2*N0DT0T+1 
ELSE IF(J.LE.LP2) THEN 
L=LIN(J) 

ELSE 

L=LIN(J) +NODTOT 
END IF 

DO 70 I=K, IA-2 

Ad, J)=F06(LL, L) 

LL=LL+1 

70 CONTINUE 

60 CONTINUE 

LL=1 

DO 80 I=K, IA-2 

B(I)=B1 (LL) 

LL=LL+1 

SO CONTINUE 

C 

C FINALLY WE WRITE THE C MATRIX 

C 

DO 90 I=I,LMAX-1 
C(I)=0. 

90 CONTINUE 

C(LMAX)=1. 

C 

C NOW WE ARE READY TO CALL THE LP SUBROUTINE 

CALL ZX3LP(A, lA. B, C, LMAX. Ml , M2, S, PSOL, DSOL, RW, IW, lER) 
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C 

c 

c 

c 



99 

C 

C 

C 



NOW TO MASSAGE OUTPUT OF ZX3LP SO THAT IT CAN BE USED 
IN THE MAIN PROGRAM OUTPUT ROUTINE 

DO 99, I=1,LMAX-1 
K=LIN<I) 

IF(I.LE.LP2) THEN 

HMNT (K) =PSOL ( I ) -HORZMO 

ELSE 

VMNT (K) =PSOL ( I ) -VERTMO 

END IF 
CONTINUE 

PP=PSOL(LMAX) 

NOW TO RETURN TO THE MAIN PROGRAM 

RETURN 

END 



SUBROUTINE EONS 



INCLUDE ’COMMON. FOR’ 

C SUBROUTINE TO CALCULATE AND STORE THE NODAL EQUATIONS OF THE GRID 

C 

C FIRST STEP. — INITIALIZE THE ARRAYS 

NODTOT=NODESH*NODESV 
DO 18 I=l,NODEEQ 

B1 ( I)=0. 

DO 17 J=1,N0DT0T*2+1 
FOGd, J)=0. 

17 CONTINUE 

18 CONTINUE 
C 

C NOW TO CALCULATE THE NODAL EQUATIONS AND WRITE THEM TO THE ROWS OF 

C 

IF(HLEN.GE. VLEN) THEN 
XMULT=HLEN 

ELSE 

XMULT=VLEN 

END IF 

M=N0DESH+2 
N=NODTOT-NODESH- 1 
L=1 



DO 19 I=M,N 

IF(JJ (I) .EQ.NODESH.OR. JJ<I> .EQ. 1) GO TO 19 
R=0. 

C CALCULATE THE COEFFICIENT OF HORZ MO. AT NODE I 

HC=2. /HLEN 
R=R+2. *HORZMO/HLEN 
FOG(L, I)=HC 
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C CALCULATE THE COEFFICIENT OF HORZ MO. TO LEFT OF NODE I 

IF(INT(HMNT(I-1) ) .EQ.O) GO TO 20 
HC=-1. /HLEN 
R=R-HORZMO/HLEN 
FOG(L, I-1)=HC 

20 CONTINUE 

C CALCULATE THE COEFFICIENT OF THE HORZ MO. TO RIGHT OF NODE I 

IF<INT<HMNT(I+1) ) .EQ.O) GO TO 21 
HC=-1. /HLEN 
R=R-HORZMO/HLEN 
FOG(L, I+1)=HC 

21 CONTINUE 

C CALCULATE THE COEFFICIENT OF THE VERT MO. AT NODE I 

VC=2. /VLEN 
R=R+2*VERTM0/VLEN 
Kl=NODTOT+I 
F0G(L,K1)=VC 

C CALCULATE THE COEFFICIENT OF THE VERT MO. ABOVE NODE I 

K2=I-N0DESH 

IF(INT<VMNT(K2) ) .EQ.O) GO TO 22 

VC=-1. /VLEN 

R=R-VERTMO/VLEN 

K3=N0DT0T+K2 

F0G(L,K3)=VC 

22 CONTINUE 

C CALCULATE THE COEFFICIENT OF THE VERT MO. BELOW NODE I 

K2=I+N0DESH 

IF(INT(VMNT(K2) ) .EQ.O) GO TO 23 

VC=-1. /VLEN 

R=R-VERTMO/VLEN 

K3=N0DT0T+K2 

F0G(L,K3)=VC 

23 CONTINUE 

K2=2*NODTOT+l 

FOG (L, K2) =-PTLOAD ( I ) 

B1 (L)=R 

C CHECK TO SEE IF RHS OF EQUATION IS NEGATIVE, BRANCH ACCORDINGLY 

C 

IF(R.LT.O. > 60 TO 24 
DO 25 J=1,2*N0DT0T+1 

FOG(L, J)=F06<L, J)*XMULT 

25 CONTINUE 

B1 <L)=B1 <L)*XMULT 
GO TO 26 

24 DO 27 J=1,N0DT0T*2+1 

F06(L, J)=F06(L, J)*-XMULT 
27 CONTINUE 

B1 (L) =B1 (L) *-XMULT 

26 CONTINUE 
L=L+1 

19 CONTINUE 

RETURN 
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END 

SUBROUTINE CONST 
INCLUDE 'COMMON. FOR' 



C 

C THIS SUBROUTINE WRITES THE ADDITIONAL INEQUALITY CONSTRAINT 

C EQUATIONS SELECTED BY THE USER TO ENSURE THAT THE BENDING 

C MOMENT IN A HORIZONTAL OR VERTICAL BEAM DOES NOT EXCEED 

C THE BEAM PLASTIC BENDING MOMENT <IE MOMENT BETWEEN NODES) 

C THIS ALLOWS THE FORMATION OF PLASTIC HINGES AT LOCATIONS OTHER 

C THAN THE NODES 

C THIS ROUTINE ALSO HAS THE CAPABILITY TO CONSTRAIN MOMENTS OF A 

C MEMBER AS HAVING ALREADY REACHED THE YIELD MOMENT. THAT IS 

C SIMILAR TO HAVING A MEMBER THAT HAS ALREADY FAILED OR HAS MINIMAL 

C STRENGTH LEFT DUE TO EXTREME LOCAL LOADING. THIS IS 

C ACCOMPLISHED BY MAKING THE ADDITIONAL CONSTRAINT EQUATION 

C AN EQUALITY CONSTRAINT RATHER THAN AN INEQUALITY CONSTRAINT 

C RE- INITIATE HMNT ?t VMNT ARRAYS FOR USE AS MULTIPLIERS TO 

QEQST=0 DETERMINE EX I STANCE OF PARTICULAR MOMENTS BASED ON BOUNDARIES 
5 CONTINUE 

10 IF(IAUT0.EQ.2) THEN 

PRINT*, ' INPUT THE REFERENCE NODE FOR THE ADDITIONAL' 

PRINT*. 'CONSTRAINT, AND A 1 FOR HORIZONTAL MEMBER' 

PRINT*. 'OR A 2 FOR VERTICAL MEMBER (INTEGERS)' 

PRINT*. 'ENTER ZEROS FOR NO CONSTRAINTS' 

READ <5,*) NN.MANDIR 
IF(NN.EQ.O) GO TO 690 

PRINT*. ' INPUT THE DISTANCE FROM THE REFERENCE NODE' 

PRINT*, 'TO LOCATE THE CONSTRAINT (REAL) ' 

READ (5,*) XX 

BHKLOC (NN)=XX 
BVKLOC(NN)=XX 

PRINT*, 'ENTER A 1 TO INDICATE THAT THE MEMBER HAS A' 

PRINT*, 'FULLY DEVELOPED YIELD MOMENT ALREADY DEVELOPED' 

PRINT*, 'NOTE: THIS TYPE MUST BE ENTERED LAST' 

PRINT*, 'ENTER A O IF THE HINGE IS NOT ALREADY FORMED' 

PRINT*, ' (INTEGERS) ' 

READ (5,*) NQ 
NEQST=NEQST+NQ 
GO TO 20 
END IF 

DO 700 NN=1, NODTOT 

20 BHPVAL=0.0 

BHPLOC=0. O 
BVPVAL=0.0 
BVPL0C=0.0 
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BHLAV=0.0 
BHLBV=0. O 
BVLAV=0.0 
BVLBV=0.0 



IF(IAUT0.EQ.2) THEN 

BHKMNT ( NN ) =HORZMO 
BVKMNT (NN) =VERTMO 

END IF 



DO 100 J=1,LPT 

IF (LDIND(J) .EQ.NN. AND. DILOC(J) 
GO TO 100 

IF (LDIND(J) .EQ.NN. AND. DILOC(J) 
.AND. BHKMNT (NN) .GE.HORZMO) 



.EQ.0.0) 



.GT. 0.0 
THEN 



BHPVAL=D1VAL(J) 
BHPL0C=D1L0C (J) 



ELSE IF (LDIND(J) .EQ.NN. AND. DILOC(J) .LT. 0.0 
+ .AND. BVKMNT (NN) .GE.VERTMO) THEN 



B VPVAL=ABS ( D 1 V AL ( J ) ) 
B VPLOC=ABS ( D 1 LOG ( J ) ) 



END IF 

100 CONTINUE 



DO 200 J=1,LNLDN 

IF(LD2ND(J) .EQ.NN .AND. ((D2BVL(J) .GT. 0.0. OR. 

+ D2AVL(J) .GT.0.0) ) .AND. BHKMNT (NN) . GE. HORZMO) 

+ THEN 

BHLAV=D2AVL(J) 

BHLBV=D2BVL(J) 

ELSE IF (LD2ND(J) .EQ.NN .AND. (D2BVL(J) .LT. O. 

+ .0R.D2AVL(J) .LT.0.0) . AND. BVKMNT (NN) 

+ .GE.VERTMO) THEN 

BVLAV=ABS (D2AVL ( J ) ) 

BVLBV=ABS (D2BVL ( J ) ) 



END IF 

200 CONTINUE 

C IF THERE ARE NO LOADS ASSOCIATED WITH THIS NODE, CONTNUE 

IF(BHPVAL.EQ.0.0 .AND. BVPVAL .EQ. 0.0 .AND. BHLAV .EQ. 
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+ 0.0 .AND. BHLBV .EQ. 0.0 .AND. BVLAV .EQ. 0.0 .AND. 

+ BVLBV .EQ. 0.0) GO TO 680 

IF (BHPVAL. EQ. O. AND. BHLAV. EQ. 0. O. AND. BHLBV. EQ. O. 0) 

+ GO TO 450 

C CHECK FOR VERTICAL OR HORIZONTAL FOR MANUAL MODE 
IF<MANDIR.EQ.2) GO TO 450 
C HORIZONTAL BEAM CALCULATIONS 
XL=HLEN 
A=BHPLOC 
Q=BHPVAL 

STEP=HLEN/ 1 000 . O 
XX=BHKLOC(NN) 



C INCREASING OR CONSTANT LINE LOAD WITH MAX MOMENT TO THE LEFT OF 
C OR AT A POINT LOAD 

IF (BHLAV. LE. BHLBV. AND. XX.LE. A) THEN 
P=BHLBV-BHLAV 
W=BHLAV 



X1=(1.-XX/XL)*HMNT (NN) 

X2= (XX/XL) *HMNT (NN+1 ) 

X3= ( 1 . -A/XL) *Q*XX+. 5*W*XL*XX+P*XL*XX/6. O 
+ -.5*W*XX**2-P*XX**3/ (6.0*XL) 

X4=H0RZM0* ( 1 . + ( 1 . -XX/XL) *HMNT (NN) + (XX/XL) 

+ *HMNT(NN+1> ) 

IDIR=1 

XDCNST (NN)=XX 
IF (NQ. EQ. 1) THEN 
NDCNST(NN) =-NN 
ELSE 

NDCNST(NN)=NN 
END IF 

CALL C0NST2 ( XI , X2, X3, X4, IDIR. NN) 

C INCREASING OR CONSTANT LINE LOAD WITH MAX MOMENT 
C TO RIGHT OF POINT LOAD OR NO POINT LOAD 

ELSE IF (BHLAV. LE. BHLBV. AND. XX. GT. A) THEN 

P=BHLBV-BHLAV 

W=BHLAV 

XXX=XL-XX 



Xl= (XXX/XL) *HMNT (NN) 

X2= ( 1 . -XXX/XL) *HMNT (NN+1 ) 

X3=(A/XL) *Q*XXX+.5*W*XL*XXX+P*XL*XXX/3. 
+-. 5*W*XXX**2 -.5*P*XXX**2+P*XXX**3/ (6.0*XL) 
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X4=H0RZM0*( 1 . + ( 1 . -XXX/XL) ♦HMNT(NN+1 ) -*-(XXX/XL) 
♦HMNT (NN) ) 

IDIR=1 

XD IF(NQ.EQ. DTHEN 

NDCNST(NN)=-NN 
ELSE 

NDCNST (NN) =NN 
ENDIF 



CALL C0NST2 ( X 1 , X2, X3, X4, IDIR, NN) 



C DECREASING LINE LOAD WITH MAX MOMENT TO THE LEFT OF OR AT 
C A POINT LOAD OR WITH NO POINT LOAD 

ELSE IF (BHLBV.LT.BHLAV. AND. XX.LE. A) THEN 
P=BHLAV-BHLBV 
W=BHLBV 

X1=(1.-XX/XL)*HMNT(NN) 

X2= (XX/XL) *HMNT (NN+1 ) 

X3= ( 1 . -A/XL) *Q*XX+. 5*W*XL*XX+P*XL*XX/3. 5*W*XX**2 
+ +P*XX**3/ (6. 0*XL) 

X4=H0RZM0*- ( 1 . + ( 1 . -XX/XL) *HMNT (NN) + ( XX/XL) *HMNT(NN+1 ) ) 
IDIR=1 

XDCNST (NN) =XX 
IF(NQ.EQ. DTHEN 
NDCNST (NN)=-NN 
ELSE 

NDCNST (NN)=NN 
ENDIF 

CALL CONST2 ( X 1 , X2, X3, X4 , IDIR, NN) 

C DECREASING LINE LOAD WITH MAX MOMENT TO THE RIGHT OF THE 
C POINT LOAD 

ELSE IF(BHLBV.LT.BHLAV. AND. XX.GT. A) THEN 

P=BHLAV-BHVBV 

W=BHLBV 

XXX=XL-X 

Xl= (XXX/XL) *HMNT(NN) 

X2=( 1 . -XXX/XL) *HMNT (NN+1 ) 

X3=(A/XL) *Q*XXX+. 5*W*XL*XXX+P*XL*XXX/6. -. 5*W+XXX**2 
+ +P*XXX»*3/ (6.0+XL) 

X4=H0RZM0* ( 1 . + ( 1 . -XXX/XL) *HMNT (NN+1 ) + (XXX/XL) +HMNT (NN) ) 
IDIR=1 

XDCNST (NN)=XX 
IF(NQ.EQ. DTHEN 
NDCNST (NN)=-NN 
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ELSE 

NDCNST (NN)=NN 
END IF 

CALL C0NST2 (X 1 , X2, X3, X4, IDIR, NN) 

END IF 

C VERTICAL MEMBER CALCULATIONS 
450 CONTINUE 

C SKIP THIS SECTION IF THERE ARE NO LOADS ON THE VERT MEMBER 

IF ( (BVPLOC. EQ. O. O. OR. BVPVAL. EQ. O. 0) . AND. BVLAV. EQ. 0. O 
+ . AND.BVLBV.EQ.0.0) GO TO 680 

C SKIP THIS SECTION IN MANUAL IF HORIZONTAL CONSTRAINT 

IF(MANDIR.EQ. 1) GO TO 680 
XL=VLEN 
A=BVPLOC 
Q=BVPVAL 
XX=BVKLOC(NN) 



C INCREASING LINE LOAD WITH MAX MOMENT TO THE LEFT OF 
C OR AT A POINT LOAD 

IF (BVLAV. LE.BVLBV. AND. XX.LE. A) THEN 
P=BVLBV-BVLAV 
W=BVLAV 

X1=(1.-XX/XL)*VMNT(NN) 

X2= ( X X / XL ) *VMNT ( NN+NODESH > 

X3= ( 1 . -A/XL) *Q*XX+. 5*W*XL*XX+P*XL*XX/6. -. 5*W*XX**2 
+ -P*XX**3/ (6. 0*XL) 

X 4= VERTMO* <1. + (1.-XX/XL) * VMNT ( NN ) + ( X X / XL ) * VMNT ( NN+NODESH ) ) 
IDIR=2 

XDCNST (NN)=-XX 
IF(NO.EQ. DTHEN 
NDCNST (NN)=-NN 
ELSE 

NDCNST (NN)=NN 
END IF 

CALL C0NST2 ( X 1 , X2, X3, X4, IDIR, NN) 

C INCREASING LINE LOAD WITH MAX MOMENT TO RIGHT OF 
C POINT LOAD OR NO PT LOAD 

ELSE IF(BVLAV.LE.BVLBV. AND. XX.GT. A) THEN 

P=BVLBV-BVLAV 

W=BVLAV 

XXX=XL-XX 
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X1=(XXX/XL) *VMNT (NN) 

X2= ( 1 . -XXX/XL) *VMNT (NN+NODESH) 

X3= (A/XL) *Q*XXX+. 5*W»XL*XXX+P*XL*XXX/3. -.5*W*XXX**2 
-. 5*P»XXX**2+P*XXX**3/ (6. 0*XL) 

X4=VERTM0* ( 1 . + ( 1 . -XXX/XL) *VMNT (NN+NODESH) 

+(XXX/XL) *VMNT(NN) ) 

IDIR=2 



XDCNST(NN)=-XX 
IF (NQ.EQ. 1 ) THEN 
NDCNST(NN)=-NN 
ELSE 

NDCNST(NN)=NN 
END IF 



CALL C0NST2 (X 1 , X2, X3, X4, IDIR, NN) 



C DECREASING LINE LOAD WITH MAX MOMENT TO THE LEFT OF OR AT 
C A POINT LOAD OR WITH NO POINT LOAD 

ELSE IF (BVLBV.LT.BVLAV. AND. XX.LE. A) THEN 
P=BVLAV-BVLBV 
W=BVLBV 

XL) *VMNT (NN+NODESH) 

X3= ( 1 . -A/XL) *Q*XX+. 5*W*XL*XX+P*XL*XX/3. -. 5*W*XX*+2 
+ +P*XX**3/ (6.0*XL) 

X4=VERTM0* (l. + (l.-XX/XL) *VMNT (NN ) + ( X X / XL ) *VMNT ( NN+NODESH ) ) 
IDIR=2 

XDCNST(NN)=-XX 
IF (NO. ED. 1) THEN 
NDCNST(NN) =-NN 
ELSE 

NDCNST(NN)=NN 
END IF 

CALL C0NST2 ( X 1 , X2, X3, X4, IDIR, NN) 

C DECREASING LINE LOAD WITH MAX MOMENT TO THE RIGHT OF THE 
C POINT LOAD 

ELSE IF(BVLBV.LT.BVLAV. AND. XX.GT. A) THEN 

P=BVLAV-BVVBV 

W=BVLBV 

XXX=XL-X 

Xl= (XXX/XL) *VMNT(NN) 

X2= ( 1 . -XXX/XL) *VMNT (NN+NODESH) 

X3= (A/XL) *Q*XXX+. 5*W*XL*XXX+P*XL*XXX/6. -. 5*W*XXX**2 
+ +P*XXX**3/ (6.0*XL) 

X4=VERTM0* ( 1 . + ( 1 . -XXX/XL) *VMNT (NN+NODESH) 
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680 



690 



+ + (XXX/XL) *V«NT (NN) ) 

IDIR=2 

XDCNST(NN)=-XX 
IF(NQ.EQ. DTHEN 
NDCNST(NN)=-NN 
ELSE 

NDCNST(NN)=NN 
END IF 

CALL C0NST2 (XI , X2, X3, X4, IDIR, NN) 

END IF 

I F ( BHPVAL . EQ . 0 . O . AND . BVPV AL . EQ . O . O . AND . 

+ BHLAV. EQ. O. 0. AND. BHLBV. EQ. O. 0. AND. 

BVLAV. EQ. O. O. AND. BVLBV. EQ. 0.0. AND. 

+ IAUT0.EQ.2) THEN 

IF(MANDIR.EQ. DTHEN 
XMULT1=HMNT(NN) 

XMULT2=HMNT ( NN+ D 
XL=HLEN 
XMO=HORZMO 
IDIR=1 

ELSE IF(MANDIR.EQ.2)THEN 
XMULT1=VMNT (NN) 

XMULT2=HMNT (NN+NODESH) 

XL=VLEN 

XMO=VERTMO 

IDIR=2 

END IF 

X1=(1.-XX/XL)*XMULT1 
X2=(XX/XL) *XMULT2 
X 3=0.0 

X4=XM0* ( 1 . + ( 1 . -XX/XL) *XMULT1 + XX/XL*XMULT2) 
IFdDIR.EQ. DTHEN 
XDCNST(NN)=XX 
ELSE IFdDIR.EQ. 2) THEN 
XDCNST(NN)=-XX 
END IF 

IF(NQ.EQ.l) THEN 

NDCNST (NN)=-NN 

ELSE 

NDCNST (NN)=NN 
END IF 

CALL C0NST2(X1, X2, X3, X4, IDIR.NN) 

END IF 



IF(IAUT0.EQ.2) THEN 
PRINT*, ' 

PRINT*, 'ENTER 1 FOR ADDITIONAL CONSTRAINTS OR 0' 



700 

720 

500 

10 

15 

20 

S— 

30 

35 

40 
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PRINT*, 'FOR NO ADDITIONAL RESTRAINTS REQUIRED' 
READ (5,*) nil 
IF<IIII.EQ. 1) 60 TO 10 
GOTO 720 
END IF 

CONTINUE 

RETURN 

END 

SUBROUTINE C0NST2 ( X 1 , X2, X3, X4, ID IR, NDREF) 

INCLUDE 'COMMON. FOR' 

L=2*N0DT0T 

IEQNM=IEQNM+1 

N2=NDREF+1 

IF(IDIR.EQ.2) GO TO 20 
DO 10 I=1,NDEF-1 

A2(IEQNM, I)=0. 

CONTINUE 

A2 ( lEQNM, NDREF) =X1 
A2<IEQNM,N2)=X2 
DO 15 I=N2+1,L 

A2(IEQNM, I)=0. 

CONTINUE 

A2( IEQNM,L+1 >=X3 
B2<IEQNM)=X4 

60 TO 40 

CONTINUE 
J=NODTOT+NDREF 
K=J+NODESH 
DO 25 1=1, J-1 

A2<IEQNM, I)=0. 

CONTINUE 

A2(IEQNM, J)=X1 
DO 30 I=J+1,K-1 

A2<IEQNM, I)=0. 

CONTINUE 

A2<IEQNM,K)=X2 
DO 35 I=K+1,L 

A2(IEQNM, I)=0. 

CONTINUE 

A2(IEQNM.L-*-l)=X3 

B2(IEQNM)=X4 

CONTINUE 

RETURN 

END 
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SUBROUTINE BEAMCK 
INCLUDE 'COMMON. FOR' 

CHARACTER* 9 NAME . NAME 1 , NAME2 , NAME3 , NAME4 , NAMES 
CHARACTER* 9 NAME6, NAME7, NAMES 

NAME1=' I /NODE#' 

NAME='LD2ND #' 

NAME2='D2AVL' 

AME5=' 

NAME6='BVLAV' 

NAME7='BVLBV' 

NAME8=' NNN/LNLDN' 

DO 999 I=1,LNLDN 
999 CONTINUE 



DO 700 1=1, NODTOT 

BHKMNT<I)=0.0 

BHKL0C(I)=0.0 

BVKMNT(I)=0.0 

BVKL0C(I)=0.0 



BHPVAL=0. O 
BHPL0C=0.0 
BVPVAL=0.0 
BVPL0C=0.0 



BHLAV=0. O 
BHLBV=0. O 
BVLAV=0.0 
BVLBV=0.0 



5000 



FORMAT (2X,A9, IX, 14, 5X , A9, G12. 2, 3X , A9, G12. 
G12.2,3X, A9,G12.2) 



2,3X, A9, 



DO 100 J=1,LPT 

IF (LDIND(J) .EQ. I. AND. DILOC(J) .EQ.0.0) GO TO lOO 
IF (LDIND(J) .EQ. I. AND. DILOC(J) .GT. 0.0) THEN 

BHP VAL=D 1 VAL ( J ) *PP 
BHPL0C=D1L0C(J) 



ELSE IF (LD1ND<J) .EQ. I. AND.DILOC(J) .LT.0.0) THEN 
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BVPVAL=ABS < D 1 VAL ( J ) ) *PP 
BVPLOC=ABS (DILOC ( J ) ) 



ENDIF 



100 



CONTINUE 



DO 200 NNN=1,LNLDN 

IF( (LD2ND<NNN) .EQ. I) .AND. (D2BVL(NNN) .GT. 0.0. OR. 

+ D2AVL(NNN) .GT.0.0) ) THEN 

BHLAV=D2AVL (NNN) *PP 
BHLBV=D2BVL (NNN) *PP 

ELSE IF ( (LD2ND (NNN) . EQ. I ) . AND. (D2BVL (NNN) . LT. O. O. OR. 
+ D2AVL(NNN) .LT.0.0) ) THEN 

BVLAV=ABS (D2AVL (NNN) ) *PP 
BVLBV=ABS ( D2BVL ( NNN ) ) *PP 

ENDIF 

200 CONTINUE 

C IF THERE ARE NO LOADS ASSOCIATED WITH THIS NODE, CONTNUE 

IF(BHPVAL.EQ.0.0 .AND. BVPVAL .EQ. 0.0 .AND. BHLAV .EQ. 

+ 0.0 .AND. BHLBV .EQ. 0.0 .AND. BVLAV .EQ. 0.0 .AND. 

+ BVLBV .EQ. 0.0) GO TO 450 



IF (BHPVAL. EQ. O. AND. BHLAV. EQ. 0. O. AND. BHLBV. EQ. O. O) 
+ GO TO 450 

C HORIZONTAL BEAM CALCULATIONS 
XL=HLEN 
A=BHPLOC 
Q=BHPVAL 
BM1=HMNT(I) 

BM2=HMNT(I+1) 

STEP=HLEN / 1 OOO . O 

BMXA=0. O 
BMXB=0. O 
BMX=0.0 

DO 300 NNN= 1,1000 
X=NNN*STEP 
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C 

C 



C 



300 

450 

C 

C 



XX = X 

SET X FOR SECOND INTERVAL IF TO RI6HT OF POINT LOAD 
IF(A.NE.O.O. AND.Q.NE.0.0. AND.X.6T. A> THEN 
X=X-A 
END IF 

INCREASINB OR CONSTANT LINE LOAD 
IF(BHLAV.LE.BHLBV) THEN 
P=BHLBV-BHLAV 
W=BHLAV 

RC0M=(BM2-BM1 ) *X/XL+W*XL*X*. 5+ 

+ P*XL»X/6.0-.5*W*X*»2-P*X**3/ 

+ (6.0*XL) 

DECREASING LINE LOAD 
ELSE 

P=BHLAV-BHLBV 

W=BHLBV 

RCOM= <BM2-BM1 ) *X/XL+W*XL*X*. 5+ 

+ P*XL*X/3. 0-. 5*W*X**2-P*X**2*. 5 

+P*X**3/ (6. 0*XL) 

END IF 



IF(Q.EQ.O.O.OR. A.EQ.0.0) THEN 
BMX=RC0M+Bt11 
ELSE IF (X.LE.A) THEN 

BMXA=RCOM+Bm+ ( 1-A/XL) *Q*X 

ELSE 

BMXB=BMXA+RCOM-A*Q*X/XL 

END IF 

IF<ABS(BMX) .LT.ABS(BMXA) .OR. ABS(BMX) 
+ .LT. ABS(BMXB) ) THEN 

IF(ABS<BMXA) .6T. ABS(BMXB) ) THEN 
BMX=BMXA 
ELSE 

BMX=BriXB 
END IF 
END IF 

IF(ABS(BHKMNT<I) ) .LT. ABS(BMX) ) THEN 
BHKMNT(I)=BMX 
BHKLOC(I)=XX 
END IF 

CONTINUE 



CONTINUE 

SET MAX MOMENT TO CONSTRAINED VALUE IF MEMBER HAS EXTRA 
CONSTRAINT EQUATION FOR IT AND CONSTRAINT IS LIMITING 



IF<XDCNST(I) .GT.0.0. AND. <BHKMNT(I) .GE.HORZMO 



n n 
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+ .OR.NDCNST(I) .LT.O) )THEN 

BHKMNT ( I ) =HORZMO 
BHKLOC ( I ) =XDCNST ( I ) 

END IF 

VERTICAL MEMBER CALCULATIONS 

SKIP THIS SECTION IF THERE ARE NO LOADS ON THE VERT MEMBER 

IF( (BVPLOC.EQ.O.OR.BVPVAL.EQ.O.O) -AND. (BVLAV. EQ. 0. 0 
+ .AND.BVLBV.EQ.0.0) ) SO TO 610 



A=BVPLOC 
Q=BVPVAL 
BM1=VMNT<I> 
BM2=VMNT ( I+NODESH) 
STEP= VLEN / 1 000 . O 
XL=VLEN 



BMXA=0.0 
BMXB=0. 0 
BMX=0.0 



DO 600 NNN= 1,1000. 

X=NNN*STEP 
XX = X 

C SET X FOR SECOND INTERVAL IF TO RIGHT OF POINT LOAD 
IF<A.NE.0.0. AND.Q.NE.0.0. AND. X.6T. A> THEN 



C INCREASING OR CONSTANT LINE LOAD 
IF (BVLAV. LE.BVLBV) THEN 
P=BVLBV-BVLAV 
W=BVLAV 



C 



+ 

+ 



DECREASING 

ELSE 



END IF 



RCOM= (BM2-BM1 ) «■ ( X/XL> +W*XL*X*. 5+ 
(P*XL*X/6. O) - ( . 5*W*X**2) - 
(P*X*-»-3/ <6.0*XL) ) 

LINE LOAD 

P=BVLAV-BVLBV 

W=BVLBV 

RC0M=(BM2-BM1) *X/XL+W*XL*X*.5+ 

P*XL*X/3. 0-. 5*W*X**2-P*X**2*. 5 
+P*X**3/ (6. 0*XL) 



IF(Q.EQ.O.O.OR. A.EQ.0.0) THEN 
BMX=RC0M+BM1 
ELSE IF (X.LE.A) THEN 



n n 
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600 

610 



700 



C 

C 



BMXA=RC0M+BM1+ ( 1-A/XL) *Q*X 

ELSE 

BMXB=BMXA+RCOM-A*Q*X/XL 

END IF 

IF(ABS(BMX) .LT. ABS(BMXA) .OR. ABS(BMX) 

+ -LT. ABS(BMXB) ) THEN 

IF(ABS(BMXA) .GT. ABS(BMXB) ) THEN 
BMX=BliXA 
ELSE 

BMX=BMXB 
ENDIF 
END IF 

IF(ABS(BVKMNT<I) ) .LT. ABS(BMX) ) THEN 
BVKMNT ( I ) =BMX 
BVKLOC ( I ) =X 
ENDIF 

CONTINUE 

CONTINUE 

IF(XDCNSTd) .LT.0.0 .AND. (BVKMNT(I) . GE. VERTMO 
+ -OR. NDCNST(I) .LT.0.0) ) THEN 

BVKMNT ( I ) =VERTf10 
BVKLOC < I ) =ABS ( XDCNST < I ) ) 

ENDIF 

CONTINUE 

RETURN 

END 



SUBROUTINE OUTPUT 

INCLUDE ’COMMON. FOR’ 

CHARACTER* 15 SIDE, BBCC , D I R , FLAG 
CHARACTER*20 CNST 

THIS ROUTINE OUTPUTS THE GRID LOADING SCHEME AND ALL MOMENTS AND LO. 
FACTORS 

CHECK ITERATION NUMBER IF IT IS ’1’ THEN PRINT GRID DATA 
IF ( ITERAT.NE. 1 ) GO TO 50 
PRINT OUT GRID SIZE 
WRITE (15,1000) 

WRITE (15,1100) NODESH.HLEN 
WRITE (15,1 200 ) NODES V , VLEN 



C 



PRINT OUT THE GRID BOUNDARY CONDITIONS 
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WRITE (15,1300) 

DO HEN 

IBOUND=ITOP 

SIDE='TOP’ 

ELSE IF (I.EQ.2) THEN 
IBOUND=IBTM 
S I DE=' BOTTOM’ 

ELSE IF (I.EQ.3) THEN 
IBOUND=ILFT 
SIDE=’LEFT’ 

ELSE 

IB0UND=IR6T 

SIDE=’RI6HT’ 

END IF 

IF (IBOUND .EQ. 1) THEN 

BBCC=’ SIMPLE SUPORTED’ 

ELSE 

BBCC=’ CLAMPED’ 

END IF 

WRITE (15,1400) SIDE.BBCC 

15 CONTINUE 

PRINT*, HORZMO, VERTMO 
C WRITE THE ZX3LP ERROR NUMBER 

WRITE(15, 1425) HORZMO, VERTMO 

1425 FORMAT(////, ’ HORZ BEAM PLASTIC BENDING MOMENT = ’,G12.2, 
+//,’ VERT BEAM PLASTIC BENDING MOMENT = ’,G12.2) 



WRITE (15, 1450) lER 

1450 FORMAT (//, 5X, ’ ZX3LP ERROR NUMBER = ’,14) 

C PRINT OUT ALL THE POINT LOADS ON THE GRID 

WRITE (15,1500) 

WRITE (15,1600) 

DO 25 1=1, LPT 

IF(ABS(D1L0C(I) ) -LT.O.Ol) THEN 
DIR=’ON THE NODE’ 

ELSE IF (DILOC DIR=’TO THE RIGHT’ 

ELSE 

DIR=’BELOW’ 

END IF 

WRITE (15,1700) LDIND(I) ,DIR,ABS(D1L0C(I) ) ,D1VAL(I) , I 
25 CONTINUE 



C PRINT OUT LINE LOADS 

WRITE (15,1800) 
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DO 30 I=1,LNLDN 
IF (D2BVL(I) .GT. 0-0) THEN 
NB=LD2ND ( I)+l 

ELSE 

NB=LD2ND ( I ) +NODESH 

ENDIF 

WRITE (15,1900) LD2ND<I) ,NB,ABS(D2AVL<I) ) ,ABS(D2BVL(I) ) , I 
30 CONTINUE 



1000 FORMAT ( ' 1 ' , 5X , ' GRID CHARACTERISTICS AND LOADING',///) 
1100 FORMAT (6X, 'NUMBER OF NODES HORIZONTAL = ’,I3,5X, 
+'HORIZONTAL LENGTH =',F6.2) 

1200 FORMAT (6X, 'NUMBER OF NODES VERTICAL = ',I3,5X, 

+'VERTICAL LENGTH =',F6.2) 

1300 FORMAT (' O' , 5X ,' BOUNDARY COND ITIONS' ,/ ) 

1400 FORMAT (6X , A6, T25, A15) 

1500 FORMAT (///, 6X ,' POINT LOADS ADDED TO THE GRID',/) 

1600 FORMAT (3X,'REF NODE' , T14, ' DIR FM N0DE',T35, 

+' DISTANCE FROM NODE' , T57, ' MAGNITUDE' , T70, 

+' POINT LOAD NUMBER') 

1700 FORMAT (6X, 13, T14, A13, T45, F6. 2, T5S, FIO. 2, T72, 13) 

1800 FORMAT (8 < / ) , 5X , ' LINE LOADS' ,//, 6X, ' BETWEEN NODES', 

+T30, 'REF NODE MAG' , T45, ' OTHER NODE MAG', 

+T65, 'LINE LOAD NUMBER',/) 

1900 FORMAT (6X, 13, ' AND ' , 13, T28, FIO. 2, T44, FIO. 2, T70, 13) 



40 CONTINUE 

C PRINT OUT THE FINAL NODAL LOADING THAT CALCULATIONS ARE BASED ON 

WRITE (15,2300) 

2300 FORMAT (8 ( / ) , 5X, ' NODE LOADS USED IN CALCULATIONS',//) 

WRITE (15,2400) 

2400 FORMAT (6X,'NODE NUMBER' , T20, ' NODE LOAD’) 

DO 50 I=1,N0DT0T 

WRITE ( 15,2500) I , PTLOAD ( I ) 

2500 FORMAT ( lOX, 13, T21 , F8. 2) 

50 CONTINUE 

WRITE (15, 2550) ITERAT 

C PRINT OUT THE EXTRA CONSTRAINT DATA 

WRITE (15,2525) 

2525 FORMAT (//,2X, 'ADDITIONAL CONSTRAINTS USED’,/) 

IFdAUTO.EQ. 1) THEN 

WRITE (15, 2527) 

ELSE IF( IAUT0.EQ.2)THEN 
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2527 

2528 

2529 

2540 



57 

2545 

C 

2550 

2600 

2700 

2800 

60 

2900 



WRITE (15, 2528) 

ELSE 

WRITE (15, 2529) 

ENDIF 

FORMAT (2X, 'AUTOMATIC CONSTRAINTS WERE USED') 

FORMAT (2X, 'MANUAL CONSTRAINTS WERE USED') 

FORMAT (2X, 'NO ADDITIONAL CONSTRAINTS WERE USED') 

WRITE (15, 2540) 

FORMAT ( 6X , ' BETWEEN NODES ' , 3X , ' LOCAT ION', 

TYPE USED' ) 

DO 57, I=1,N0DT0T 

N0D1=ABS (NDCNST ( I ) ) 

IF(XDCNSTd) .GT.0.0) THEN 
N0D2=ABS (NDCNST ( I ) ) +1 
ELSE IF(XDCNSTd) .LT.0.0) THEN 
N0D2=ABS (NDCNST ( I ) ) +NODESH 
ELSE 

60 TO 57 
ENDIF 

IF(NDCNSTd) .6T.0.0) THEN 
CNST=' INEQUALITY ' 

ELSE 

CNST=' EQUALITY ' 

ENDIF 

WR ITE( 15, 2545) NODI, N0D2,ABS(XDCNST(I) ) , CNST 
CONTINUE 

FORMAT (6X, 13, ' AND ' , 13, 6X , 610. 2, A20) 

PRINT OUT MOMENTS AT EACH NODE AND THE LOAD FACTOR 

FORMAT (2X,//, '***** ITERATION NUMBER IS' , 14, ' *****',//) 
WRITE (15,2600) 

FORMAT (2 (/) ,6X, 'MOMENTS AND LOAD FACTOR') 

WRITE (15,2700) 

FORMAT (//,6X,'N0DE NUMBER' , T25, ' HOR I Z MNT' , T40, ' VERT MNT' ) 
DO 60 I=1,N0DT0T 

WRITE (15, 2800) I , HMNT ( I ) , VMNT ( I ) 

FORMAT (8X, 13, T26, E8. 2, T42, E8. 2) 

CONTINUE 

WRITE (15,2900) PP 

FORMAT (///,6X,'THE GRID LOAD FACTOR = ',68.2) 



C 



PRINT OUT MAX MOMENT AND LOCATION FOR EACH MEMBER 



3000 

4050 

4100 

70 

4300 

80 

4200 

5000 
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WRITE (15, 3000) 

FORMAT (3 (/), ////,6X, 'MAXIMUM MOMENTS AND LOCATION 
+ ON LOADED MEMBERS’,/) 

WRITE (15, 4050) 

FORMAT (6X, 'BETWEEN NODES’ , T30, ’MA(3NITUDE’ ,T45, 

+ 'LOCATION') 

WRITE(15,4100) 

FORMAT (6X, 'HORIZONTAL MEMBERS’ ) 

DO 70 I=1,N0DT0T 

IF(BHKMNTd) .EQ.0.0) 60 TO 70 
IF(ABS(BHKMNT(I) ) .GT. (HORZMO+ . 01 ) ) THEN 
FLAG=’**’ 

ELSE 

FLAG=' ' 

ENDIF 

WRITE (15, 4200) I , I +1 , BHKMNT ( I ) , BHKLOC ( I ) , FLAG 
CONTINUE 



WRITE (15, 4300) 

FORMAT (6X, 'VERTICAL MEMBERS') 

DO 80 I=1,N0DT0T 

IF(BVKMNTd) .EQ.0.0) GO TO 80 
IF(ABS(BVKMNT( I ) ) .GT. (VERTMO+.Ol) ) THEN 
FLAG =’♦*’ 

ELSE 

FLAG =' 

ENDIF 

WR I TE ( 1 5 , 4200 ) I , I +NODESH , B VKMNT (I ) , B VKLOC ( I ) , FL AG 
CONTINUE 
WRITE (15,5000) 

FORMAT (6X, 13, ’ AND ' , 13, T28, 610. 4, 8X , GIO. 4, 2X , A15) 
FORMAT (////, 5X, INDICATES > YIELD MOMENT') 

RETURN 

END 



- 95 - 



APPENDIX B 



IMSL LINEAR PROGRAMMING ROUTINE 
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The linear progr ammi ng routine used by GRIDS is taken 
from the International Mathematics and Statistics Library (IMSL). 
6RIDS uses routine ZX3LP- The attached pages describe the 
routine input requirements and output formatting and variables- 
This information is provided for a better understandi ng of the 
dimensioning process that GRIDS performs. These pages have been 
reproduced from the IMSL library reference manual . 
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PURPOSE 


- Solve the linear programmi ng problem via the revise 
simplex algorithm — easy to use version 


USAGE 


“• CALL ZX3LP (A, IA,B,C,N,M1,M2,S,PS0L,DS0L,RW, IW, lER) 



ARGUMENTS 



A 


- Matrix of dimension Ml+M2-»-2 by N containing the 
coefficients ot the Ml rows followed by the 
coefficients of the M2 equality constraints (input) 
The last two rows are used only as working storage. 


lA 


- Row dimension of matrix A exactly as specified 
in the dimension statement in the calling 
program, (input) Two rows of A are required for 
working storage, and therefore, lA must not be 
less than M1+M2+2. 


B 


Vector of length Ml+M2->*2 containing the right 
hand sides of the inequality constrai nts. ( i nput ) 
The last two elements of B are used as working 
storage. 


C 


Vector of length N containing the coefficients 
of the objective function, (input) 


N 


h4umber of unknowns in the model- (input) 


Ml 


Number of inequality constraints, (input) 


M2 


Number of equality constraints- (input) 


S 


Value of the objective function- (output) 


PSOL - 


Vector of length N containing the primal 
solution, (output) PSOL is also used as work 
storage and therefore must have length at least 
MAX (n,Ml+M2) . 


DSOL - 


Vector of length M1+M2+2 containing the DUAL 
solution, (output) That is, DSOL ( 1 ) , 

DSOL (M1+M2) contain the solution to the problem 
MIN BT*Y subject to AT*Y is greater than or 
equal to C and Y greater than or equal to 0 
where AT = A-transpose and BT = B— transpose- 
When the primal problem has equality 
constraints, the correspondi ng components of the 
dual solution are unconstrained. DSOL (M1+M2+1 ) 

and DSOL (M1+M2+2) are used as working storage- 


RW 


Work Vector of length (M1+M2+2) * (M1+M2+2) + 3*M1 
2*M2+4. 



-98- 



IW - Work Vector of length 2*M2 3*M1 4. 



lER - ERROR indicator, (output) 

Terminal Error 

lER = 130 Indicates that lA is less than 
M1H-M2+2. 

lER = 131 Indicates that the cost criterion 
has unbounded values. 

lER = 132 Indicated that the Maximum number 

of iterations was reached in ZXOLP. 

lER = 133 Indicated that no feasible solution 
ex i sts. 

Warnig (with fix) 

lER = 70 Indicates that some artifical 

variables remain in the solution 
basis at a zero level after phase 
1. This condition can be caused by 
having redundant constraints. 
Nevertheless, a solution is computed 
and returned in PSOL and DSOL. 

Algor i thm 



To solve the linear programming problem. 



Maximize C^PSOL^ + ... C^PSOL^^ = SP 

subject to 





a . , PSOL , 
1 1 1 


... 




a . . .PSOL. , 
1 N N 


< 


8. i — 1,..., 

1 




a. .PSOL, 
1 1 1 


... 




a... PSOL., 
iN N 


= 


B. i=M +1,. 

1 1 










PSOL . 
J 


- 


0 j > • - • ? ^ 


where M 


= Ml + M2. 












The DUAL 


1 i near progr ammi ng 


problem is, 








Mi ni mi ze 


bjDSOLj 




• • • ^ ®M 


DSOL,,, = SD 
N 


subject 


to 














a, DSOL, 
Ij 1 


+ ... 






> 


C. j = l 1 










DSOL. 

1 


> 


■ 

■ 

■ 

II 

•l-l 

o 


DSOL. 

1 


unrtestr i cted in sign 


when 




i = Mj + 1 M 
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ZX3LP computes the solution to the primal problem, PSOL, the 
solution to the dual problem, DSOL, and the values oi the 
objective -function S=SP=SD. 

ZX3LP calls ZXOLP which solve the linear programming problem 
by the revised simplex method- 



where M = Ml M2. 
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APPENDIX C 

SAMPLE PROBLEMS AND RESULTS 
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EXAMPLE 1 



GRID CHARACTERISTICS AND LOADING 



NUMBER OF NODES HORIZONTAL = 6 HORIZONTAL LENGTH = 10.00 

NUMBER OF NODES VERTICAL = 5 VERTICAL LENGTH = 10.00 

BOUNDARY CONDITIONS 

TOP CLAMPED 

BOTTOM CLAMPED 

LEFT CLAMPED 

RIGHT CLAMPED 



HORZ BEAM PLASTIC BENDING MOMENT = O. lOE+05 

VERT BEAM PLASTIC BENDING MOMENT = 0. lOE+05 



ZX3LF ERROR NUMBER = O 



POINT LOADS ADDED TO THE GRID 

REF NODE DIR FM NODE DISTANCE FROM NODE MAGNITUDE 
8 ON THE NODE 0.00 1.00 



LOAD NUMBER 
1 



LINE LOADS 



BETWEEN NODES 



REF NODE MAG 



OTHER NODE MAG 



LINE LOAD NUMBER 
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NDDE LOADS USED IN CALCULATIONS 



NODE NUMBER 
1 



4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 

22 

OT 

24 

25 

26 

27 

28 

29 

30 



NODE LOAD 
O. 00 
0.00 
0.00 
O. 00 
0.00 
0.00 
0.00 
1.00 
O. 00 
0.00 
0.00 
0 . 00 
O. 00 
0.00 
0.00 
0 . 00 
O. 00 
0.00 
O. 00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 



ITERATION NUMBER IS 1 ***** 



ADDITIONAL CONSTRAINTS USED 

NO ADDITIONAL CONSTRAINTS WERE USED 

BETWEEN NODES LOCATION TYPE USED 



MOMENTS AND LOAD FACTOR 
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NODE NUMBER 

1 

O 

X— 

T 

♦_/ 

4 

5 

6 

7 

8 
7 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 

23 

24 

25 

26 

27 

28 

29 

30 



HORIZ MNT 
O. OOE+00 
O. OOE+OO 
O. OOE+OO 
O. OOE+OO 
O. OOE+OO 
O. OOE+OO 
-. lOE+05 
O. 10E+05 
-. lOE+05 
-. 40E+04 
-.70E+04 
-. 10E+05 
-. lOE+05 
0.84E+04 
0.68E+04 
-. lOE+05 
--89E+04 
-. lOE+05 
-. lOE+05 
-. lOE+05 
-. lOE+05 
-.55E+04 
-. lOE+05 
-. lOE+05 
O. OOE+OO 
O. OOE+OO 
0. OOE+OO 
O. OOE+OO 
O. OOE+OO 
O. OOE+OO 



THE GRID LOAD FACTOR = O. 80E+04 



MAXIMUM MOMENTS AND LOCATION ON 

BETWEEN NODES MAGNITUDE 

HORIZONTAL MEMBERS 
VERTICAL MEMBERS 



VERT MNT 
0. OOE+OO 
-. 10E+O5 
-. lOE+05 
-. lOE+05 
-. lOE+05 
0. OOE+OO 
O. OOE+OO 
0. lOE+05 
0.30E+04 
-. lOE+05 
-. lOE+05 
O. OOE+OO 
0. OOE+OO 
-. 10E+05 
-. lOE+05 
-. lOE+04 
-. lOE+05 
O. OOE+OO 
O. OOE+OO 
-. lOE+05 
-. 7SE+04 
lOE+05 
-.78E+04 
0. OOE+OO 
O. OOE+OO 
-. 10E+05 
-. lOE+05 
-. lOE+05 
-. lOE+05 
O. OOE+OO 



LOADED MEMBERS 
LOCATION 



** INDICATES > YIELD MOMENT 



BBHBi 

i|!|j 




:S|Pf 


1 I'M. 


d in 


lit' 








i-i* 

li 


.isf 4 -- m"t\ 

.joioioiot m - raap^'lilM'i 

Ksyr-- 


U|,;r 

w 




*1- i 


l^ifc 


W”s 


Pj 


0^^-' 




EJ 
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EXAMPLE 2 



GRID CHARACTERISTICS AND LOADING 



NUMBER OF NODES HORIZONTAL = 5 HORIZONTAL LENGTH = 10.00 

NUMBER OF NODES VERTICAL = 5 VERTICAL LENGTH = 10.00 

BOUNDARY CONDITIONS 



TOP 

BOTTOM 

LEFT 

RIGHT 



SIMPLE SUPORTED 
SIMPLE SUPORTED 
SIMPLE SUPORTED 
SIMPLE SUPORTED 



HORZ BEAM PLASTIC BENDING MOMENT = O. lOE+05 

VERT BEAM PLASTIC BENDING MOMENT = O. lOE+05 



ZX3LP ERROR NUMBER = O 



POINT LOADS ADDED TO THE GRID 

REF NODE DIR FM NODE DISTANCE FROM NODE MAGNITUDE LOAD NUMBER 



LINE LOADS 



BETWEEN 


NODES 


REF NODE MAG 


OTHER NODE MAG 


LINE LOAD NUMBER 


7 AND 


8 


1.00 


1.00 


1 


7 AND 


12 


1 . 00 


1 . 00 




3 AND 


13 


1.00 


1.00 


“T 


12 AND 


13 


1.00 


1.00 


4 



-105- 



NODE LOADS USED IN CALCULATIONS 



NODE 



NUMBER 


NODE LOAD 


1 


0.00 


o 

X. 


0-00 


3 


0-00 


4 


0-00 


5 


0- 00 


6 


0. 00 


7 


10.00 


8 


10.00 


9 


0-00 


10 


0- 00 


11 


0. 00 


12 


10.00 


13 


10. 00 


14 


0.00 


15 


0. 00 


16 


0.00 


17 


0. 00 


18 


0.00 


19 


0- 00 


20 


0 . 00 


21 


0.00 


22 


0.00 




0. 00 


24 


0.00 


25 


0.00 



**** ITERATION NUMBER IS 1 ***** 



ADDITIONAL CONSTRAINTS USED 

AUTOMATIC CONSTRAINTS WERE USED 

BETWEEN NODES LOCATION TYPE USED 



MOMENTS AND LOAD FACTOR 



NODE NUMBER 
1 

3 

4 

5 

6 



HORIZ MNT 
O. OOE+00 
O. OOE+00 
O. OOE+00 
0. OOE+00 
O. OOE+00 
0. OOE+00 



VERT MNT 
O.OOE+00 
O. OOE+00 
O.OOE+00 
O. OOE+00 
O. OOE+OO 
O.OOE+00 
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7 


0-89E+04 


0. lOE+05 


8 


0. lOE+05 


0. lOE+05 


9 


0-33E+04 


0-67E+04 


10 


0.00E^“00 


0- OOE-^00 


11 


0. OOE+00 


0. OOE+00 


12 


0- lOE+05 


0- lOE+05 


13 


0- lOE+05 


0. lOE+05 


14 


0-33E+04 


0- lOE+05 


15 


0- OOE+00 


0. OOE+00 


16 


0. OOE+00 


0. OOE+00 


17 


0.78E+04 


0.22E+04 


18 


0. lOE+05 


llE+04 


19 


0.98E-03 


0. lOE+05 


20 


0- OOE+00 


0- OOE+00 


21 


0- OOE+00 


0. OOE+00 


22 


0. OOE+00 


0. OOE+00 


O-T 


0. OOE+00 


0. OOE+00 


24 


0. OOE+00 


0- OOE+00 


25 


0. OOE+00 


0. OOE+00 



THE GRID LOAD FACTOR = O. 18E+03 



MAXIMUM MOMENTS AND 


LOCATION ON 


LOADED MEMBERS 




BETWEEN NODES 


MAGNITUDE 


LOCATION 




HORIZONTAL MEMBERS 


7 AND 8 


0. 1170E+05 


5.620 




12 AND 13 


0. 1222E+05 


5.000 




VERTICAL MEMBERS 


7 AND 12 


0. 1222E+05 


5.000 




8 AND 13 


0. 1222E+05 


5. 000 





** INDICATES > YIELD MOMENT 



***♦ ITERATION NUMBER IS 2 ***** 



ADDITIONAL CONSTRAINTS USED 
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AUTOMATIC CONSTRAINTS WERE USED 



BETWEEN NODES 


LOCATION 


TYPE USED 


7 AND 12 


5. 0 


INEQUALITY 


8 AND 13 


5-0 


INEQUALITY 


12 AND 13 


5-0 


INEQUALITY 


MOMENTS AND 


LOAD FACTOR 




NODE NUMBER 


HORIZ MNT 


VERT NNT 


1 


0- OOE+00 


0- OOE+00 


2 


0. OOE+00 


0- 00E-»“O0 




0- OOE+00 


0- OOE+00 


4 


0- OOE+00 


0- OOE+00 


5 


0- OOE+00 


0- OOE+00 


6 


0- OOE+00 


0- OOE+00 


7 


0-79E+04 


0-62E+04 


8 


0-66E+04 


0-89E+04 


9 


-- 17E+04 


0- lOE+05 


10 


O.OOE+00 


0. OOE+00 


11 


0- OOE+00 


0- OOE+00 


12 


0-61E+04 


0-99E+04 


13 


0. lOE+05 


0-71E+04 


14 


0-50E+04 


0- lOE+05 


15 


0- OOE+00 


0- OOE+00 


16 


0- OOE+00 


0- OOE+00 


17 


0- lOE+05 


-.68E-4-02 


18 


0- lOE+05 


-- 14E+04 


19 


--29E-02 


0. lOE+05 


20 


0. OOE+00 


0. OOE+00 


21 


O.OOE+00 


0- OOE+00 




0- OOE+00 


0- OOE+00 


O-T 


0. OOE+00 


O.OOE+00 


24 


0- OOE-4-00 


0- OOE+00 




0- OOE+00 


0. OOE+00 



THE GRID LOAD FACTOR = 0. 16E+03 



MAXIMUM MOMENTS AND LOCATION ON LOADED MEMBERS 

BETWEEN NODES MAGNITUDE LOCATION 

HORIZONTAL MEMBERS 



M Nj 
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7 AND 


8 


0. 1058E+05 


2.910 


12 AND 


13 


0. lOOOE+05 


5. 000 


VERTICAL 


MEMBERS 






7 AND 


12 


0. lOOOE+05 


5 . 000 


8 AND 


13 


0. lOOOE+05 


5.000 



INDICATES > YIELD MOMENT 
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EXAMPLE 3 



GRID CHARACTERISTICS AND LOADING 



NUMBER OF NODES HORIZONTAL = 5 HORIZONTAL LENGTH = 10.00 

NUMBER OF NODES VERTICAL = 5 VERTICAL LENGTH = 10.00 

BOUNDARY CONDITIONS 

TOP CLAMPED 

BOTTOM CLAMPED 

LEFT CLAMPED 

RIGHT CLAMPED 



HORZ BEAM PLASTIC BENDING MOMENT = 0. lOE+05 

VERT BEAM PLASTIC BENDING MOMENT = O. lOE+05 



ZX3LP ERROR NUMBER = O 



POINT LOADS ADDED TO THE GRID 

REF NODE DIR FM NODE DISTANCE FROM NODE MAGNITUDE LOAD NUMBER 
7 ON THE NODE 0.00 1.00 1 



LINE LOADS 



BETWEEN NODES 



REF NODE MAG 



OTHER NODE MAG 



LINE LOAD NUMBER 
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NODE LOADS USED IN CALCULATIONS 



NODE NUMBER 
1 

2 

"T 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 
16 

17 

18 

19 

20 
21 

22 

o*r 

X— 

24 

25 



NODE LOAD 
0 . 00 
0.00 
0.00 
0.00 
0.00 
0.00 
1 . OO 
0.00 
0 . 00 
0.00 
0.00 
0.00 
0.00 

0.00 
0.00 
O. 00 
0.00 
0.00 
O. 00 
0.00 
O. 00 
0.00 
O. 00 



ITERATION NUMBER IS 1 ***** 



ADDITIONAL CONSTRAINTS USED 

NO ADDITIONAL CONSTRAINTS WERE USED 

BETWEEN NODES LOCATION TYPE USED 



MOMENTS AND LOAD FACTOR 



NODE NUMBER 

1 

'-i 

T 

4 

5 

6 



HORIZ MNT 
O.OOE+00 
O. OOE+00 
O- OOE+00 
O-OOE+00 
O.OOE+00 
lOE+05 



VERT MNT 
O. OOE+00 
lOE+05 
-- lOE+05 
lOE+05 
O. OOE-^00 
O. OOE+00 
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7 


0. lOE+05 


0. lOE+05 


8 


-. lOE+05 


0.80E+04 


9 


-. lOE+05 


-. lOE+05 


10 


-. lOE+05 


0. OOE+OO 


1 1 


-. lOE+05 


0. OOE+OO 


12 


0.20E+04 


-. lOE+05 


13 


-. lOE+05 


0.60E+04 


14 


-.80E+04 


-. lOE+05 


15 


-. lOE+05 


0. OOE+OO 


16 


-. lOE+05 


0. OOE+OO 


17 


-. lOE+05 


-. 60E+O4 


18 


-. 20E+04 


-. lOE+05 


19 


-. lOE+05 


-.60E+04 


20 


-. lOE+05 


0. OOE+OO 


21 


0. OOE+00 


0. OOE+OO 




0. OOE+00 


-. lOE+05 


OT 


0. OOE+00 


-. lOE+05 


24 


0. OOE+OO 


-. lOE+05 


25 


0. OOE+00 


0. OOE+OO 



THE GRID LOAD FACTOR = 0.80E+04 



MAXIMUM MOMENTS AND LOCATION ON LOADED MEMBERS 

BETWEEN NODES MAGNITUDE LOCATION 

HORIZONTAL MEMBERS 
VERTICAL MEMBERS 



** INDICATES > YIELD MOMENT 
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EXAMPLE 4 



GRID CHARACTERISTICS AND LOADING 



NUMBER OF NODES HORIZONTAL = 4 HORIZONTAL LENGTH = 10.00 

NUMBER OF NODES VERTICAL = 4 VERTICAL LENGTH = 10.00 

BOUNDARY CONDITIONS 



TOP 

BOTTOM 

LEFT 

RIGHT 



SIMPLE SUPORTED 
SIMPLE SUPORTED 
SIMPLE SUPORTED 
SIMPLE SUPORTED 



HORZ BEAM PLASTIC BENDING MOMENT = O. lOE+05 

VERT BEAM PLASTIC BENDING MOMENT = O. 10E+05 



ZX3LP ERROR NUMBER = O 



POINT LOADS ADDED TO THE GRID 

REF NODE DIR FM NODE DISTANCE FROM NODE MAGNITUDE LOAD NUMBER 



LINE LOADS 

BETWEEN NODES 

2 AND 6 

h AND 10 

5 AND 6 

6 AND 7 



REF NODE MAG 

1.00 
1 . 00 
1 . 00 
1 . 00 



OTHER NODE MAG 

1.00 

1.00 

1.00 

1 . 00 



LINE LOAD NUMBER 

1 

2 

3 

4 
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NODE LOADS USED IN CALCULATIONS 



NUMBER 


NODE LOAD 


1 


0- 00 


o 


5.00 


3 


0 . oo 


4 


0. 00 


5 


5. 00 


6 


20. 00 


7 


5. OO 


8 


0.00 


7 


0. 00 


10 


5. 00 


11 


0.00 


12 


0- 00 


13 


0. 00 


14 


0. 00 


15 


0.00 


16 


0 . oo 



ITERATION NUMBER IS 1 ***** 



ADDITIONAL CONSTRAINTS USED 

AUTOMATIC CONSTRAINTS WERE USED 

BETWEEN NODES LOCATION TYPE USED 



MOMENTS AND LOAD FACTOR 



NODE NUMBER 


HORIZ MNT 


VERT MNT 


1 


0. OOE+00 


0. OOE+OO 


O 

J— 


0. OOE+00 


0. OOE+OO 


3 


O.OOE+00 


0. OOE+OO 


4 


0. OOE+00 


0. OOE+OO 


5 


0. OOE+OO 


0 . OOE+OO 


6 


0. lOE+05 


0. 10E+05 


~T 

J 


-. 50E+03 


0. lOE+05 


8 


0. OOE+OO 


0. OOE+OO 


9 


0. OOE+OO 


0. OOE+OO 


10 


0. 10E+05 


0. 45E+04 
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11 


0. lOE+05 


-. 98E-03 


12 


0. OOE+00 


0. OOE+00 


13 


O.OOE+00 


O.OOE+00 


14 


0. OOE+00 


0. OOE+00 


15 


O.OOE+00 


0. OOE+OO 


16 


O.OOE+00 


0. OOE+OO 



THE GRID LOAD FACTOR = O. 18E+03 



MAXIMUM MOMENTS AND LOCATION ON LOADED MEMBERS 



BETWEEN NODES 
HORIZONTAL MEMBERS 

5 AND 6 

6 AND 7 

VERTICAL MEMBERS 

2 AND 6 

6 AND 10 



MAGNITUDE 

0. lOOOE+05 
9998. 

O. lOOOE+05 
0. 1034E+05 



LOCATION 

10 . 00 

O. lOOOE-01 



10 . 00 

1.940 



** INDICATES > YIELD MOMENT 



**** ITERATION NUMBER IS 2 



ADDITIONAL CONSTRAINTS USED 
AUTOMATIC CONSTRAINTS WERE USED 



BETWEEN 


NODES 


LOCATION 


TYPE USED 


5 AND 


6 


10. 


INEQUALITY 


6 AND 


10 


1.9 


INEQUALITY 



MOMENTS AND 


LOAD FACTOR 




NODE NUMBER 


HORIZ MNT 


VERT MNT 


1 


0. OOE+OO 


O. OOE+OO 


2 


0. OOE+OO 


0. OOE+OO 






-115- 



3 


O.OOE+OO 


O.OOE+OO 


4 


O.OOE+00 


O.OOE+OO 


5 


0. OOE+OO 


O.OOE+OO 


6 


0. 10E+05 


0. lOE+05 


7 


0. 13E+04 


0. lOE+05 


a 


0. OOE+00 


O.OOE+OO 


9 


O.OOE+OO 


O.OOE+OO 


lO 


0. lOE+05 


0.27E+04 


11 


0.65E+04 


0.35E+04 


12 


0. OOE+OO 


O. OOE+OO 


13 


0. OOE+OO 


0. OOE+OO 


14 


O. OOE+OO 


0. OOE+OO 


15 


0. OOE+OO 


0. OOE+OO 


16 


0. OOE+OO 


0. OOE+OO 



THE GRID LOAD FACTOR = 0. 18E+03 



MAXIMUM MOMENTS AND LOCATION ON LOADED MEMBERS 



BETWEEN NODES 
HORIZONTAL MEMBERS 

5 AND h 

6 AND 7 

VERTICAL MEMBERS 

2 AND 6 

6 AND 10 



MAGNITUDE 



O. lOOOE+05 
O. lOOOE+05 

O. lOOOE+05 
O. lOOOE+05 



LOCATION 

10 . 00 

O. 1400 

10 . 00 
1 . 940 






** INDICATES > YIELD MOMENT 
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EXAMPLE 5 

GRID CHARACTERISTICS AND LOADING 



NUMBER OF NODES HORIZONTAL = 3 HORIZONTAL LENGTH = 10.00 

NUMBER OF NODES VERTICAL = 3 VERTICAL LENGTH = 10.00 

BOUNDARY CONDITIONS 



TOP 

BOTTOM 

LEFT 

RIGHT 



SIMPLE SUPORTED 
SIMPLE SUPORTED 
SIMPLE SUPORTED 
SIMPLE SUPORTED 



HORZ 


BEAM 


PLASTIC 


BENDING 


MOMENT = 


30 


VERT 


BEAM 


PLASTIC 


BENDING 


MOMENT = 


30 



ZX3LP ERROR NUMBER = 0 

POINT LOADS ADDED TO THE GRID 

REF NODE DIR FM NODE DISTANCE FROM NODE MAGNITUDE LOAD NUMBER 



LINE LOADS 
BETWEEN NODES 
4 AND 5 



REF NODE MAG 
1.00 



OTHER NODE MAG 

1.00 



LINE LOAD NUMBER 
1 
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NODE LOADS USED IN CALCULATIONS 



NODE NUMBER 
1 
X. 

3 

4 

5 

6 

7 

8 
9 



NODE LOAD 
0.00 
0.00 
0.00 
5.00 
5. 00 
0.00 
O. 00 
0.00 
0.00 



**** ITERATION NUMBER IS 1 ***** 



ADDITIONAL CONSTRAINTS USED 

AUTOMATIC CONSTRAINTS WERE USED 

BETWEEN NODES LOCATION TYPE USED 



MOMENTS AND LOAD FACTOR 



NODE NUMBER 
1 



4 

5 

6 

7 

8 
9 



HORIZ MNT 
O.OOE+00 
O. OOE+00 
O. OOE+00 
O. OOE+00 
0.30E+02 
0. OOE+00 
0. OOE+00 
0. OOE+00 
O. OOE+00 



VERT MNT 
0. OOE+00 
O.OOE+00 
O. OOE+OO 
0. OOE+OO 
O. 30E+02 
O. OOE+OO 
0. OOE+OO 
0. OOE+OO 
O. OOE+OO 



THE GRID LOAD FACTOR = 2.4 



MAXIMUM MOMENTS AND LOCATION 



ON LOADED MEMBERS 




i 



-1 



BETWEEN NODES MAGNITUDE LOCATION 

HORIZONTAL MEMBERS 

4 AND 5 46.87 6.250 

VERTICAL MEMBERS 



** INDICATES > YIELD MOMENT 



**** ITERATION NUMBER IS 2 ***** 



ADDITIONAL CONSTRAINTS USED 

AUTOMATIC CONSTRAINTS WERE USED 

BETWEEN NODES LOCATION TYPE USED 
4 AND 5 6.3 INEQUALITY 



MOMENTS AND LOAD FACTOR 



NODE NUMBER 

1 

T> 

3 

4 

5 

6 

7 

8 
9 



H0RI2 MNT 
O. OOE+00 
0. OOE+00 
O. OOE+00 
O.OOE+00 
O. 15E+02 
O. OOE+00 
O. OOE+00 
O. OOE+00 
O. OOE+00 



VERT MNT 
O.OOE+00 
O. OOE+00 
O. OOE+00 
O. OOE+00 
0. 30E+02 
O. OOE+OO 
O. OOE+OO 
O. OOE+OO 
O. OOE+OO 



THE GRID LOAD FACTOR = 1.8 



MAXIMUM MOMENTS AND LOCATION ON LOADED MEMBERS 



BETWEEN NODES MAGNITUDE 

HORIZONTAL MEMBERS 

4 AND 5 30.00 

VERTICAL MEMBERS 



LOCATION 

6.250 



** 



** INDICATES > YIELD MOMENT 
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